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Abstract 

Monkeypox  virus  (MPV)  is  a  zoonotic  Orthopoxvirus  and  a  potential  biothreat  agent  that  causes  human  disease 
with  varying  morbidity  and  mortality.  Members  of  the  Orthopoxvirus  genus  have  been  shown  to  suppress  antiviral 
cell  defenses,  exploit  host  cell  machinery,  and  delay  infection-induced  cell  death.  However,  a  comprehensive  study 
of  all  host  genes  and  virus-targeted  host  networks  during  infection  is  lacking.  To  better  understand  viral  strategies 
adopted  in  manipulating  routine  host  biology  on  global  scale,  we  investigated  the  effect  of  MPV  infection  on 
Macaca  mulatta  kidney  epithelial  cells  (MK2)  using  GeneChip  rhesus  macaque  genome  microarrays.  Functional  ana¬ 
lysis  of  genes  differentially  expressed  at  3  and  7  hours  post  infection  showed  distinctive  regulation  of  canonical 
pathways  and  networks.  While  the  majority  of  modulated  histone-encoding  genes  exhibited  sharp  copy  number 
increases,  many  of  its  transcription  regulators  were  substantially  suppressed;  suggesting  involvement  of  unknown 
viral  factors  in  host  histone  expression.  In  agreement  with  known  viral  dependence  on  actin  in  motility,  egress,  and 
infection  of  adjacent  cells,  our  results  showed  extensive  regulation  of  genes  usually  involved  in  controlling  actin 
expression  dynamics.  Similarly,  a  substantial  ratio  of  genes  contributing  to  cell  cycle  checkpoints  exhibited  con¬ 
certed  regulation  that  favors  cell  cycle  progression  in  G1,  S,  G2  phases,  but  arrest  cells  in  G2  phase  and  inhibits 
entry  into  mitosis.  Moreover,  the  data  showed  that  large  number  of  infection-regulated  genes  is  involved  in  mole¬ 
cular  mechanisms  characteristic  of  cancer  canonical  pathways.  Interestingly,  ten  ion  channels  and  transporters 
showed  progressive  suppression  during  the  course  of  infection.  Although  the  outcome  of  this  unusual  channel 
expression  on  cell  osmotic  homeostasis  remains  unknown,  instability  of  cell  osmotic  balance  and  membrane 
potential  has  been  implicated  in  intracellular  pathogens  egress.  Our  results  highlight  the  role  of  histones,  actin,  cell 
cycle  regulators,  and  ion  channels  in  MPV  infection,  and  propose  these  host  functions  as  attractive  research  focal 
points  in  identifying  novel  drug  intervention  sites. 


Introduction 

Monkeypox  vims  is  a  double-stranded  DNA  virus  and 
one  of  the  human  pathogenic  orthopoxviruses  that 
include  Variola  (VARV),  cowpox  (CPX),  and  Vaccinia 
(VACV)  viruses.  The  virus  causes  a  disease  that  mani¬ 
fests  similarly  to  smallpox,  but  with  milder  morbidity 
and  lower  mortality  rates  [1],  Variation  in  MPV  viru¬ 
lence  has  been  observed  and  mapped  to  defined  geo¬ 
graphic  origins,  e.g.,  virus  isolates  from  Central  Africa 
are  more  virulent  than  those  from  Western  Africa 
[2,3] 
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Recent  advances  in  molecular  biology  and  genomics 
have  improved  our  understanding  of  viral  infection 
and  replication  mechanisms.  Monkeypox  virus  has  a 
relatively  large  genome  of  about  196,858  base  pairs, 
encoding  190  open  reading  frames,  which  constitute 
the  bulk  of  the  material  needed  for  viral  replication  in 
cell  cytoplasm  [4].  Viral  entry  into  cells  is  dependent 
on  cell  type  and  viral  strain,  and  occur  after  an  initial 
attachment  to  cell  surface  through  interactions 
between  multiple  viral  ligands  and  cell  surface  recep¬ 
tors  [5]  such  as  chondroitin  sulfate  [6]  or  heparan  sul¬ 
fate  [7,8].  Subsequent  crossing  of  cell  membrane  is 
mediated  by  a  viral  fusion  event  with  cell  membrane 
under  neutral  pH  conditions  [9],  or  by  endosomal 
uptake  via  a  macropinocytosis-like  mechanism  that 
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involves  actin  [10,11]  and  low  pH-dependent  steps 
[12].  Once  in  the  cell  cytoplasm,  the  virus  releases  pre¬ 
packaged  viral  proteins  and  enzymatic  factors  that  dis¬ 
able  cell  defenses  and  stimulate  expression  of  early 
genes  [13-15].  Synthesis  of  early  proteins  promotes 
further  uncoating,  DNA  replication,  and  production  of 
intermediate  transcription  factors.  In  following  stage, 
intermediate  genes  are  transcribed  and  translated  to 
induce  the  expression  of  late  genes  that  function 
mainly  as  structural  proteins,  enzymes,  and  early  tran¬ 
scription  factors.  Eventually,  membrane  structures  will 
appear  and  unit  virion  genomes  processed  from  DNA 
concatemers  are  assembled  into  nascent  virions  that 
contain  all  enzymes,  factors,  and  genetic  information 
needed  for  a  new  infectious  cycle. 

The  detailed  available  information  about  viral  gene 
functions  and  its  programmed  expression  during  infec¬ 
tion  exceeds  current  knowledge  of  corresponding  events 
in  the  host.  Furthermore,  although  poxviruses  are  con¬ 
sidered  one  of  the  most  self-sufficient  viral  families,  they 
remain  unable  to  reproduce  in  extracellular  environment 
and  known  to  have  limited  host  range,  which  suggest 
dependence  on  host  elements  [16,17].  Therefore,  identi¬ 
fication  of  these  specific  host  elements  and  pathways 
that  are  essential  for  viral  replication  will  enrich  our 
knowledge  of  host  response  to  viral  infection,  and  may 
prove  valuable  in  identifying  potential  targets  for  anti¬ 
viral  therapies. 

Microarrays  have  been  used  in  genome  exploration 
and  profiling  with  special  focus  on  understanding 
dynamics  of  viral  gene  expression  and  pathogenesis 
[18,19].  However,  a  paucity  of  work  employed  this  tool 
in  examining  host  response  to  infections  with  poxviruses 
generally  [20-22],  and  more  specifically  in  the  case  of 
MPV.  Because  combining  microarray  technology  with 
modern  data  mining  tools  allows  further  information 
extraction  at  genome-wide  levels,  we  used  whole  gen¬ 
ome  rhesus  macaque  microarrays  in  combination  with 
Ingenuity  Pathways  Analysis  (Ingenuity®  Systems,  http:// 
www.ingenuity.com)  to  investigate  the  effect  of  MPV 
infection  on  host  Maccaca  mulata  kidney  epithelial  cells 
transcriptome,  and  address  gaps  in  host  response  during 
MPV  infection.  Functional  and  canonical  pathway  analy¬ 
sis  of  differentially  expressed  genes  at  3  and  7  hours 
post-infection  (hpi)  time  points  validated  many  of  the 
known  host  gene  responses  to  poxvirus  infection  and 
introduced  new  sets  of  interesting  functions  and  path¬ 
ways  in  areas  of  cell  death  and  apoptosis,  actin 
dynamics,  ion  channels  and  transport,  and  cell  cycle  reg¬ 
ulation.  Our  data  points  to  a  vital  role  for  these  cell 
functions  in  MPV  infection,  and  hence,  signify  their 
value  in  poxviruses’  infection  diagnosis  and  treatment 
studies. 


Materials  and  methods 

Cell  culture  and  viral  infection 

Monkeypox  virus-Katako  Kombe  strain  (MPV-KK)  was 
propagated  in  Vero  E6  cells  maintained  in  Eagle’s  Mini¬ 
mum  Essential  Medium  with  non-essential  amino  acids 
(EMEM/NEAA)  supplemented  with  2  mM  L-glutamine, 
10%  heat-inactivated  fetal  bovine  serum  (FBS),  10  mg/L 
Gentamycin,  250  pg/L  Fungizone,  and  buffered  at  pH 
7.4  with  10  mM  HEPES  [23].  Viral  titers  were  deter¬ 
mined  by  the  plaque  assay  [24].  As  described  previously 
[25],  monolayer  of  Vero  E6  were  inoculated  with  serial 
dilutions  of  viral  suspension  and  allowed  to  adsorb  for 
30  min  at  room  temperature.  The  viral  inoculate  was 
removed  and  cells  were  washed  twice  with  PBS  then 
incubated  with  culture  medium  for  5  days.  To  count 
plaques,  culture  medium  was  removed  and  cells  were 
fixed  and  stained  simultaneously  using  30%  formalin,  5% 
ethanol  (vol/vol)  solution  containing  1.3  g/L  crystal 
violet. 

For  time  points  and  control  samples,  MI<2  cells  were 
grown  in  the  same  culture  medium  described  above  but 
in  absence  of  antibiotics  and  Fungizon  for  at  least 
2  days  before  infection  with  MPV.  Culture  medium  was 
removed  and  cells  were  inoculated  with  crude  monkey¬ 
pox  virus-Katako  Kombe  strain  (MPV-KK)  at  MOI  of  3. 
Virus  was  adsorbed  for  30  min  at  37°C  with  gentle  rock¬ 
ing  for  15  sec  each  10  min,  then  cells  were  washed  twice 
with  room  temperature  equilibrated  PBS,  fresh  culture 
medium  was  added,  and  cells  were  incubated  for  3  or  7 
hours.  Control  cells  were  handled  identically  except  for 
exposure  to  virus. 

RNA  and  cDNA  preparation,  labeling,  hybridization,  and 
scanning 

RNA  was  extracted  using  TRIzol  LS  (Invitrogen)  accord¬ 
ing  to  the  manufacturer’s  recommended  protocol. 
Briefly,  200  pi  of  chloroform  (Sigma  Aldrich,  St.  Louis, 
MI)  was  added  to  5xl06  infected  and  non-infected  con¬ 
trol  cells  harvested  in  1  ml  of  TRIzol  LS  reagent  (Invi¬ 
trogen).  Solutions  were  mixed  thoroughly  and  incubated 
for  10-15  min  at  4°C,  then  centrifuged  for  20  min  at 
14,000  rpm  in  4°C.  After  centrifugation,  the  aqueous 
phase  was  removed,  added  to  an  equal  volume  of  2-pro- 
panol  (Sigma-Aldrich).  After  overnight  incubation  at 
-20°C,  the  mixture  was  centrifuged  for  20  min  at  14,000 
rpm  in  4°C.  Supernatant  was  removed  and  the  RNA  pel¬ 
let  was  washed  with  80%  ethanol  (Sigma  Aldrich),  dried, 
and  dissolved  in  100  pi  of  RNase-free  water  (Ambion, 
Austin,  TX).  Possibly  present  contaminating  DNA  was 
eliminated  using  Turbo  DNA-free  kit  (Ambion),  and 
RNA  clean  up  was  performed  using  RNeasy  kit  (Qiagen, 
Valencia,  CA).  Quality  and  quantity  of  RNA  was  evalu¬ 
ated  using  Experion  automated  electrophoresis  station 
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and  Experion  RNA  StdSens  analysis  kit  (Bio-Rad,  Her¬ 
cules,  CA).  Clean  sharp  peaks  representing  intact  rRNA 
were  confirmed  for  each  preparation  by  two  indepen¬ 
dent  workers. 

cDNA  was  synthesized  using  One-cycle  cDNA 
Synthesis  kit  (Affymetrix,  Santa  Clara,  CA)  in  presence 
of  poly-A  RNA  controls.  Double-stranded  cDNA  sam¬ 
ples  were  cleaned  up  and  biotin  labeling  of  antisense 
cRNA  was  carried  out  with  the  IVT  Labeling  Kit 
(Affymetrix,  P/N  901229).  Material  was  cleaned  up 
then  fragmented  before  hybridization  overnight. 
Microarrays  (Affymetrix,  P/N  900656)  were  scanned 
using  the  Affymetrix  GeneChip  scanner  following 
standard  Affymetrix  protocols  [26]  after  carrying  out 
all  washing  and  staining  steps  as  recommended  by  the 
manufacturer. 

Microarray  validation  by  RT  PCR 

RT  PCR  was  performed  utilizing  the  Superarray  human 
common  cytokines  panel  (PAHS-021E-4,  Superarray, 
Frederick  MD,  21703).  Reactions  were  performed 
according  to  manufacturer’s  guidelines  for  both  time 
points  in  10  pi  volumes  using  the  384  plate  format  that 
allow  gene  expression  analysis  in  quadruplicate.  Results 
showed  a  strong  correlation  of  gene  expression  levels 
with  that  obtained  using  microarray. 

Data  analysis 

The  Affymetrix  CEL  files  were  imported  into  Gene- 
Spring  GX  software  v  7.3.1  (Agilent),  which  allows 
multi-filter  comparisons  using  data  from  different 
experiments,  to  perform  the  normalization,  generation 
of  restriction  lists  and  functional  classifications  of  the 
differentially  expressed  genes. 

Normalization  was  applied  in  two  steps:  i)  “per  chip 
normalization”  by  which  each  measurement  was  divided 
by  the  50th  percentile  of  all  measurements  in  its  array; 
and  ii)  "per  gene  normalization”  by  which  each  treated 
sample  was  normalized  against  its  respective  control 
(mock-treated)  sample. 

The  expression  of  each  gene  was  reported  as  the  ratio 
of  the  value  obtained  relative  to  the  control  condition 
after  normalization  of  the  data.  Transcripts  whose  levels 
reproducibly  changed  were  identified  using  one-way 
parametric  analysis  of  variance  with  a  P-value  cutoff  of 
0.05  (false  discovery  rate  of  5%).  The  changes  in  tran¬ 
script  levels  are  expressed  as  the  fold  change  in  signal 
between  control  and  treated  samples. 

Two-dimensional  clustering  was  carried  out  based  on 
samples  and  genes  for  visualization  and  assessment  of 
reproducibility  in  the  profile  of  the  significant  genes 
across  biological  replicates. 


Results 

Dataset  overview 

To  characterize  and  measure  changes  in  transcriptome 
of  cells  infected  with  MPV,  we  infected  107  Macaca 
mulatta  kidney  epithelial  cells  (MK2)  and  incubated 
them  for  3  or  7  hpi.  Cells  were  infected  at  high  multipli¬ 
city  to  enhance  infection  probability  and  improve  signal 
to  noise  ratio.  Infected  cells  were  washed  after  virus 
adsorption  to  reduce  response  to  any  biological  factors 
in  the  seed,  and  to  avoid  an  infection  continuum  that 
may  desynchronize  virus  replication  stages.  We  selected 
the  two  time  points  not  to  exceed  7  hpi  to  avoid  possi¬ 
ble  entry  into  a  second  cycle  of  viral  replication,  and  to 
minimize  artifacts  that  might  arise  from  cell  lysis  after 
viral  maturation  and  release  of  intracellular  components 
[18].  We  mock-infected  a  similar  number  of  cells  and 
used  them  as  controls  for  each  time  point.  The  experi¬ 
ment  was  done  in  triplicate,  and  47,000  rhesus  macaque 
transcripts  potentially  present  in  purified  RNA  from 
each  sample  were  interrogated  by  52,000  probe  sets/ 
microarray. 

We  employed  RT-PCR  as  a  tool  to  validate  microar¬ 
ray  results.  Expression  levels  of  nine  genes  encoding 
eight  cytokines  and  one  housekeeping  protein  were 
assessed  and  the  resulting  fold  change  for  each  gene 
was  compared  with  corresponding  calculated  measure¬ 
ments  from  the  microarray  study.  The  RT-PCR  mea¬ 
surements  were  carried  out  in  quadruplicate,  and  cycle 
threshold  (Ct)  values  in  each  observation  were  normal¬ 
ized  to  the  average  of  five  housekeeping  genes  and  to 
the  expression  level  of  the  same  gene  in  control  mock- 
infected  cells  as  described  in  microarray  experiments. 
Despite  consistent  mild  lower  fold  change  values  in  all 
assessed  genes  by  RT-PCR,  a  direct  correlation  was 
evident  between  the  results  obtained  in  both  techni¬ 
ques,  and  15  out  of  the  18  gene  expression  measure¬ 
ments  in  3  hpi  (Fig.  1A)  and  7  hpi  (Fig.  IB)  were 
concordant.  The  moderately  staggered  prediction  in 
gene  regulation  in  the  other  three  measurements  coin¬ 
cided  with  large  Ct  values,  which  were  near  or  beyond 
RT-PCR  detection  sensitivity  limitations.  Likewise,  the 
observed  lower  and  consistent  fold  change  FC  esti¬ 
mates  from  RT-PCR  where  mainly  due  to  different 
normalization  methods  used  for  each  approach.  The 
strong  correlation  in  results  of  RT-PCR  and  microarray 
was  mutually  validating  and  instigated  subsequent  data 
analysis. 

MPV  compromises  host's  biological  activities  with  a 
dominating  global  downregulation 

Downregulation  was  reported  as  the  hallmark  of  gene 
expression  modulation  in  Vaccinia-iniected  human 
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Figure  1  Validation  of  microarray  results  by  quantitative  real¬ 
time  PCR.  Copy  number  of  eight  common  cytokines  and  a  house 
keeping  gene  (B2M)  were  calculated  based  on  RT-PCR  Ct  values. 
Fold  change  of  genes  expression  is  plotted  on  Y-axis  after 
normalization  to  mock-treated  samples.  Results  plotted  to  compare 
calculated  fold  change  in  expression  of  each  gene  at  3  hpi  (A)  or  at 
7  hpi  (B)  using  Microarray  (grey  bars)  and  RT-PCR  (black  bars). 


HeLa  cells  [20]  and  in  lymphocytes  of  Variola- infected 
cynomolgus  macaques  ( Macaca  fascicularis)  [21].  In 
this  study,  we  identified  2,702  transcripts  that  exhib¬ 
ited  statistically  significant  changes  in  copy  numbers 
with  P  values  <  0.05.  This  represents  about  5.7%  of  the 
total  interrogated  host  transcripts  on  the  microarray 
(Fig.  2A).  In  agreement  with  previously  reported 
results  [21,22],  the  majority  of  transcription  changes 
we  observed  were  downregulation.  More  than  89%  of 
the  regulated  genes  or  2,407  genes  exhibited  steady 
downregulation  in  both  3  and  7  hpi  time  points,  while 
only  295  genes  or  10.92%  of  total  regulated  gene 
showed  upregulation  under  the  same  statistical  con¬ 
straints  (Fig.  2B). 


Up-  and  down-regulated  genes  exhibit  varied  temporal 
regulation  distribution 

Because  MPV  exhibit  well-defined  temporal  gene 
expression  stages  that  were  classified  into  early,  inter¬ 
mediate,  and  late  stages,  we  anticipated  the  host 
response  to  exhibit  some  sort  of  an  analogous  pattern. 
This  may  be  represented  in  variable  temporal  expression 
of  regulated  host  genes,  where  a  subset  will  exhibit  a 
higher  copy  number  at  early  time  point  followed  by  a 
decline  at  late  time  point  as  the  inducing  viral  protein 
or  gene  decreases  or  disappears.  We  followed  peaks  of 
gene  expression  regulation  for  each  of  the  differentially 
expressed  genes  across  both  time  points  by  dividing 
gene  expression  fold  change  (FC)  at  late  time  point  (7 
hpi)  by  FC  of  the  same  gene  in  early  time  point  (3  hpi). 
Calculated  fold  change  ratios  (FCRs)  7/3  hpi  for  all 
2,407  down-regulated  genes  are  plotted  in  (Fig.  3A). 
Bars  represent  the  average  FCR  for  each  consecutive 
100  genes  after  sorting  them  according  to  their  FCR 
values  from  smallest  to  highest  for  data  clarity.  A  plot 
for  all  295  upregulated  genes  was  made  following  identi¬ 
cal  steps  (Fig.  3B),  but  bars  represent  the  average  FCR  of 
10  sequential  genes.  Genes  with  FCR  >  1  indicate  a  rela¬ 
tively  greater  regulation  in  later  stages,  while  <  1  ratio 
points  to  greater  regulation  in  early  stage.  Only  7.9%  of 
the  downregulated  genes  showed  higher  fold  change 
early  in  the  course  of  infection  (Fig.  3A),  while  51.1%  of 
the  upregulated  genes  exhibited  similar  temporal  expres¬ 
sion  trend  (Fig.  3B).  In  contrast  with  downregulated 
genes  (Fig.  3A),  peaks  of  host  gene  upregulation  were 
almost  equally  distributed  across  infection  time  points 
(Fig.  3B).  This  important  difference  in  the  distribution 
of  temporal  up  and  downregulated  genes  is  a  new  land¬ 
mark  of  global  host  gene  regulation  by  poxviruses  sec¬ 
ond  to  known  dominant  host  genes  downregulation 
[20-22].  While  the  steady  increase  in  host  gene  suppres¬ 
sion  denotes  known  continuous  viral  stifling  of  host 
functions  such  as  innate  immune  response  and  antiviral 
cell  defense  mechanisms,  the  variation  in  host  gene 
upregulation  intensity  seems  to  mark  the  stage  depen¬ 
dent  viral  gene  expression  and  need  for  specific  host 
material  or  machinery. 

Functional  gene  clusters 

We  used  Ingenuity  pathway  analysis  (IPA)  to  identify 
the  main  host  functions  and  canonical  pathways  influ¬ 
enced  by  viral  infection.  The  differentially  expressed 
genes  at  3,  7  hpi  time  points  fell  into  diverse  functional 
categories  including  enzymes,  transcription  regulators, 
kinases,  phosphotases,  peptidases,  transmembrane  recep¬ 
tors  and  G  protein  coupled  receptors,  transporters, 
translation  regulators,  and  micro  RNAs  (Table.  1).  To 
reduce  noise  in  functional  analysis,  we  applied  a  second 
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Figure  2  (A)  Gene  expression  overview  of  Macaca  mulatto  kidney  epithelial  cells  infected  with  MPV  (A)  Using  one-way  ANOVA,  a  total  of 
2702  elements  exhibited  reproducible  change  (P-value  <  0.05).  This  represented  5.7%  of  the  total  47000  interrogated  transcripts  on  the  GenChip. 
About  89%  of  the  differentially  expressed  transcripts  were  suppressed  and  less  than  1 1%  were  up-regulated.  (B)  Expression  data  for  significantly 
influenced  transcripts  were  hierarchically  clustered.  Columns  represent  triplicate  (A,  B,  C)  of  time  points  3  and  7  hpi.  Color  intensity  reflects  fold 
change  relative  to  control  (mock-transfected)  cells.  Red  and  green  indicate  up-  and  downregulation  respectively. 
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Sequential  number  of  binned  genes 
(100  gene/bin,  all  2407downregulated  genes) 


□  greater  upregulation  at  3  hpi 
■  greater  upregulation  at  7  hpi 
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Sequential  number  of  binned  genes 
(10  genes/bin,  all  295  upregulated  genes) 

Figure  3  Up-  and  down-regulated  transcripts  follow  different  distribution  patterns  during  infection  progression  Each  bar  represents 
average  fold-change  ratios  (FCRs)  of  gene  expression  at  7  hpi/3  hpi  for  either  100  consecutive  downregulated  genes  (A),  or  10  upregulated 
genes  (B),  after  sorting  all  genes  according  to  their  FCR  values  from  smallest  to  largest.  Genes  with  average  FCR  values  <  1  are  in  grey  and  those 
with  values  >  1  are  in  black.  Results  show  most  downregulated  genes  maintained  an  increasing  suppression  trend  in  the  course  of  infection  as 
more  than  92%  of  the  genes  gave  fold-change  ratio  >  1  (A).  The  295  upregulated  genes  exhibited  more  balanced  distribution  between  7  and  3 
hpi  with  about  51%  of  the  genes  having  higher  fold  change  at  early  time  point. 
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Table  1  Functional  distribution  of  differentially  expressed  genes  at  3,  7  hpi  time  points 


Gene  function  category 

Number  of  genes 

Up  regulated 

Down  regulated 

FC  <  1 .8 

3  hpi 

7  hpi 

3  hpi 

7  hpi 

3  hpi 

7  hpi 

transporter 

122 

10 

8 

43 

88 

69 

26 

transmembrane  receptor 

19 

3 

3 

5 

13 

11 

3 

translation  regulator 

16 

5 

5 

5 

9 

6 

2 

transcription  regulator 

261 

12 

12 

107 

197 

142 

52 

phosphatase 

48 

2 

1 

17 

33 

29 

14 

peptidase 

48 

4 

4 

9 

27 

35 

17 

other 

1278 

128 

104 

452 

852 

697 

322 

microRNA 

1 

- 

- 

1 

1 

0 

0 

ligand-dependent  nuclear  receptor 

7 

- 

- 

1 

4 

6 

3 

kinase 

134 

8 

4 

58 

96 

68 

34 

ion  channel 

11 

- 

- 

5 

9 

6 

2 

growth  factor 

8 

1 

1 

4 

4 

3 

3 

G-protein  coupled  receptor 

13 

1 

2 

8 

10 

3 

3 

enzyme 

308 

32 

31 

89 

199 

187 

78 

cytokine 

5 

2 

2 

1 

1 

2 

2 

data  filter  and  included  only  genes  that  exhibited  >1.8 
FC  in  addition  to  the  t-test  filter  with  P  values  <  0.05 
already  in  place.  A  total  of  1,013  and  1,720  genes  from 
3  and  7  hpi  time  points,  respectively,  met  analysis 
thresholds  and  were  used  in  a  comparative  functional 
analysis  to  identify  unique  and  common  influenced  host 
functions  for  both  time  points.  Results  were  scored 
based  on  Fisher’s  exact  test,  and  functions  with  P-values 
<  0.05  were  considered  statistically  significant.  When 
defined  the  -log  (P-value)  >1.301  as  a  significant  differ¬ 
ence  in  random  genes  clustering  probability  for  a  given 
function  between  the  two  time  points,  functions  clearly 
separated  into  either  time  point  dependent  or  indepen¬ 
dent  categories  (Fig.  4).  Functions  relating  to  cell  signal¬ 
ing,  cell  cycle,  cell  death,  transcriptional  modification, 
and  DNA  processing  were  more  significant  in  3  hpi 
time  point  (Fig.  4A),  and  while  protein  synthesis  and 
molecular  transport  functions  were  indentified  with  sig¬ 
nificant  P-values  only  at  7  hpi,  and  RNA  damage  and 
repair  was  found  to  be  unique  to  3  hpi  time  point  (Fig. 
4A).  Time-point  independent  functions  exhibited  com¬ 
parable  -log  (P-value)  in  both  time  points  and  included 
metabolism  of  essential  building  blocks  such  as  amino 
acids,  lipids,  and  carbohydrates  (Fig.  4B),  and  other 
functions  related  to  cell  morphology,  cellular  develop¬ 
ment,  small  molecule  biochemistry,  and  posttransla- 
tional  modification  (Fig.  4B). 

The  expression  of  histones  exhibited  interesting  pat¬ 
tern  with  potentially  important  implications  in  viral 
replication.  Functional  analysis  showed  steep  upregula- 
tion  of  histones  in  both  time  points  with  sharp  increase 
through  the  course  of  MPV  infection.  Except  for 
F1IST3F12A,  all  other  core  histone  genes  including 
HIST1H3I,  HIST1H1D,  HIST1H2BJ,  HIST1H2AD,  and 


HIST4H4  exhibited  upregulation  (Table.  2).  The  linker 
histone  H1F0,  was  downregulated  with  an  increasing 
trend.  Unlike  the  vast  increase  in  histones  expression, 
genes  encoding  major  transcription  regulators  of  his¬ 
tones  expression  such  as  CITED2,  NCOA3,  CREB1, 
YY1,  and  HDAC2  showed  increasing  suppression 
throughout  infection.  Similarly,  many  enzymes  control¬ 
ling  major  modifications  of  histones  and  chromatin 
organization  dynamics  of  the  host,  including  FBXOll, 
PRMT3,  MYST2,  MYCBP2,  and  RARS2  showed  steady 
downregulation  (Table.  2). 

Another  set  of  genes  encoding  mainly  transmembrane 
proteins  that  function  as  ion  channels  or  transporters 
were  progressively  suppressed  during  MPV  infection. 
Out  of  10  identified  channels  that  exhibited  statistically 
significant  downregulation  in  the  both  data  sets,  eight 
were  downregulated  at  1.8-fold  change  or  more  (Table. 
3).  The  majority  of  the  identified  channels  were  loca¬ 
lized  to  the  plasma  membrane  and  only  two  had 
unknown  locations  in  the  cell.  Identified  infection-regu¬ 
lated  channels  covered  the  transport  of  the  principal 
ions  Na+,  I<+,  CP,  and  Ca++  which  function  in  a  con¬ 
certed  manner  to  define  transmembrane  potential  and 
maintain  cellular  osmotic  balance. 

Canonical  pathways  analysis 

To  identify  the  effect  of  viral  infection  on  major  biologi¬ 
cal  processes  in  the  cell,  we  examined  potential  enrich¬ 
ment  of  all  differentially  expressed  genes  in  known 
canonical  pathways.  The  two  data  sets  were  analyzed 
using  IPA  library  of  canonical  pathways.  Genes  exhib¬ 
ited  >  1.8  FC  and  associated  with  a  canonical  pathway 
in  the  Ingenuity  pathways  knowledge  database  were 
considered  for  the  analysis.  Both  data  sets  contained 
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Figure  4  Comparative  functional  analysis  of  differentially  expressed  genes  at  3  and  7  hpi  (A)  Represents  functions  exhibiting  -log  (P- 
values)  differences  >  1.301  between  3  and  7  hpi  time  points.  Few  identified  functions  were  unique  to  one  time  point  only  under  adopted 
statistical  constraints.  (B)  Equally  present  functions  in  both  time  points  showing  a  -log  (P-value)  differences  <  1.301.  All  functions  were  identified 
with  P-values  <  0.05  or  -log  (P-value)  =  1.301  thresholds  (dotted  black  lines). 
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Table  2  List  of  differentially  expressed  genes  involved  in  histones  dynamics  in  MPV-  infected  cells 


Symbol 

Type 

Fold  Change 

3  h  7  h 

Entrez  Gene  Name 

FBX01 1 

enzyme 

-1.98 

-2.95 

F-box  protein  1 1 

PRMT3 

enzyme 

-1.81 

-3.04 

protein  arginine  methyltransferase  3 

MYST2 

enzyme 

-1.69 

-3.73 

MYST  histone  acetyltransferase  2 

MYCBP2 

enzyme 

-1.66 

-2.01 

MYC  binding  protein  2 

RARS2 

enzyme 

-1.41 

-1.61 

arginyl-tRNA  synthetase  2,  mitochondrial 

CITED2“ 

transcription  regulator 

-8.75 

-309.73 

Cbp/p300-interacting  transactivator,  with  Glu/Asp-rich  carboxy-terminal  domain,  2 

NC0A3“ 

transcription  regulator 

-2.29 

-4.76 

nuclear  receptor  coactivator  3 

CREB1“ 

transcription  regulator 

-1.91 

-2.59 

cAMP  responsive  element  binding  protein  1 

YY1“ 

transcription  regulator 

-1.58 

-2.60 

YY1  transcription  factor 

HDAC2 

transcription  regulator 

-1.38 

-1.57 

histone  deacetylase  2 

HIST3H2A 

histone 

-2.49 

-9.19 

histone  cluster  3,  H2a 

HI  FO 

histone 

-2.45 

-6.15 

HI  histone  family,  member  0 

HIST1H3I 

histone 

2.64 

22.09 

histone  cluster  1,  H3i 

HIST1H1D 

histone 

7.49 

42.62 

histone  cluster  1,  Hid 

HIST1 H2BJ 

histone 

16.57 

160.56 

histone  cluster  1,  H2bj 

HIST1  H2AD“ 

histone 

46.04 

994.02 

histone  cluster  1,  H2ad,1 

HIST4H4" 

histone 

1 964.01 

2989.9 

histone  cluster  4,  H4 

a  Average  gene  expression  values  is  presented  for  genes  with  multiple  identifiers. 


2,702  genes  from  which  2,281  genes  had  mapped  identi¬ 
ties.  The  number  of  genes  eligible  for  pathway  analysis 
was  680  and  1,134  genes  for  time  point  3  and  7  hpi, 
respectively.  The  significance  of  the  association  between 
a  subset  of  genes  in  data  sets  and  the  canonical  pathway 
was  measured  in  two  ways:  1)  ratio  of  the  number  of 
genes  from  the  data  set  that  mapped  to  the  pathway 
divided  by  the  total  number  of  genes  that  mapped  to 
the  canonical  pathway,  and  2)  Fischer’s  exact  test  P- 
value,  which  determines  the  probability  of  random 


association  between  the  genes  in  the  datasets  and  the 
canonical  pathway. 

Data  set  3  h  mapped  into  185  pathways  with  32  meet¬ 
ing  ratio  of  >  10%  and  70  meeting  a  threshold  of 
P  value  <  0.05.  Similar  analysis  of  data  set  7  h  resulted 
in  genes  mapping  to  176  pathways  with  72  achieving  a 
ratio  >  10%  and  44  having  P  values  <  0.05.  Common 
canonical  pathways  in  both  time  points  include  molecu¬ 
lar  mechanism  of  cancer,  glioma  signaling,  hepatocyte 
growth  factor  (HGF)  signaling,  cell  cycle  regulation  by 


Table  3  List  of  MK2  cell  genes  exhibited  differential  expression  upon  infection  with  MPV  and  encoded  ion  channels  or 
related  proteins 


Symbol 

Type 

Location 

Fold  Change 

Entrez  Gene  Name 

3  h 

7  h 

CLCN3  “ 

ion  channel 

Plasma  Membrane 

-2.02 

-6.68 

chloride  channel  3 

KCNMA1 

ion  channel 

Plasma  Membrane 

-2.11 

-3.94 

potassium  large  conductance  calcium-activated  channel,  subfamily  M,  alpha  member  1 

MST150 

ion  channel 

Unknown 

-2.00 

-3.63 

MSTP150 

KCMF1 

enzyme 

Unknown 

-1.77 

-2.81 

potassium  channel  modulatory 
factor  1 

SCNN1A 

ion  channel 

Plasma  Membrane 

-1.62 

-2.68 

sodium  channel,  nonvoltage-gated  1  alpha 

CUL5 

ion  channel 

Cytoplasm, 

Nucleus 

-1.69 

-2.51 

cullin  5 

KCTD20 

ion  channel 

Unknown 

-1.24 

-1.92 

potassium  channel  tetramerisation 
domain  containing  20 

TRPC1 

ion  channel 

Plasma  Membrane 

-1.87 

-1.87 

transient  receptor  potential  cation 
channel,  subfamily  C,  member  1 

CACNB3 

ion  channel 

Plasma  Membrane 

-1.31 

-1.83 

calcium  channel,  voltage-dependent, 
beta  3  subunit 

SCLT1 

transporter 

Plasma  Membrane 

-1.62 

-1.64 

sodium  channel  and  clathrin  linker  1 

TPCN1 

ion  channel 

Plasma  Membrane 

-1.34 

-1.60 

two  pore  segment  channel  1 

a  Average  gene  expression  values  are  presented  for  genes  with  multiple  identifiers. 
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B-cell  translocation  gene  (BTG)  family  of  proteins,  P53 
signaling,  cell  cycle:  G2/M  DNA  damage  checkpoint 
regulation,  cell  cycle:  Gl/S  checkpoint  regulation,  ataxia 
telangiectasia  mutated  protein  (ATM)  signaling,  ephrine 
receptor  signaling,  and  virus  entry  via  endocytic  path¬ 
ways.  Observed  order  variation  of  identified  pathways  in 
examined  time  points  when  sorted  by  their  gene  ratios 
or  P  values  underscores  the  time  dependent  viral  activity 
(Fig.  5). 

The  following  analysis  focuses  on  most  relevant  path¬ 
ways  that  were  identified  with  strong  statistical  support 
of  P-values  <  0.05  and  ratio  >  10%  in  either  time  point. 
Molecular  mechanisms  of  cancer  (MMC) 

With  370  molecules  localized  to  cell  membrane,  mito¬ 
chondrion  inner  and  outer  membrane,  cell  cytoplasm, 
and  cell  nucleus,  this  pathway  is  considered  one  of  the 
most  complex  known  pathways.  Forty-eight  genes  of 
data  set  7  hpi  clustered  to  MMC  to  cover  13%  of  total 
pathway  molecules  with  P-value  of  1.25e-5,  making  it 
the  first  pathway  influenced  by  MPV  infection.  The  3 
hpi  data  set  showed  lower  gene  ratio  of  9.5%  but  with  a 
tighter  correlation  at  lower  P-value  of  2.82e-6  (Fig.  6). 

MMC  pathway  can  be  triggered  by  diverse  stimuli 
including  hormones,  growth  factors,  cytokines,  and 
oncogenes.  Each  of  these  stimuli  is  capable  of  inducing 
a  variety  of  biological  processes  relating  to  cancer, 
including  cell  cycle  regulation  and  cell  cycle  check¬ 
points,  which  were  recognized  in  our  analysis  as  inde¬ 
pendent  pathways,  and  therefore  will  be  discussed 
separately.  The  majority  of  the  genes  specific  to  cancer 
mechanisms  pathway  are  involved  in  apoptosis  regula¬ 
tion.  Identified  modulation  in  expression  patterns  of 
these  genes  suggest  signaling  apoptosis  in  infected  cells, 
e.g.,  Noxa  or  PMAIP1  (phorbol-12-myristate-13-acetate- 
induced  protein  1)  shows  more  than  44-fold  upregula- 
tion  at  7  h  PI.  This  molecule  inactivates  the  anti-apop- 
totic  molecule  Bcl2  (B-cell  CLL/lymphoma  2)[27-29], 
which  regulates  mitochondrial  membrane  potential  by 
controlling  the  voltage  dependant  anion  channel 
(VDAC)  [30-32],  consequently  increasing  the  mitochon¬ 
drial  membrane  permeability  and  the  release  of 
cytchrome  C  that  serves  as  a  major  apoptosis  inducer 
[33].  In  addition  to  Bcl2  inactivation,  the  expression  of 
Bcl2  gene  itself  exhibited  a  significant  -2.83-fold  sup¬ 
pression,  which  would  enhance  apoptosis  signaling. 
Apoptosis-induced  cysteine-aspartic  acid  peptidases 
(caspases)  such  as  caspase  3,  which  enhance  the  degra¬ 
dation  of  Bcl2  [30]  and  trigger  subprograms  for  cell  dis¬ 
mantling  and  removal,  exhibited  3.2-fold  upregulation. 
Caspase  7,  8,  and  10  expression  remained  unchanged, 
and  only  caspase  6  [34],  which  functions  as  a  down¬ 
stream  enzyme  in  the  caspase  activation  cascade, 
showed  mild  -1.2  and  -1.62  FC  at  3,  7  hpi  time  points, 
respectively.  Other  pro-apoptotic  genes  like  BBC3  (BcL2 


binding  component  3)  or  PUMA  [35]  and  p21  protein 
(Cdc42/Rac)-activated  kinase  2  or  PAI<2  [36]  showed 
similar  suppression  trend  of  -1.85  and  -1.99  FC  at  7  hpi, 
respectively. 

Proteosomal  degradation  has  been  associated  with 
apoptosis.  Our  results  show  a  4.05  FC  upregulation  of 
nuclear  factor  of  kappa  light  polypeptide  gene  enhancer 
in  B-cells  inhibitor,  epsilon  NFKBIE,  which  mediates 
cytoplasmic  sequestering  of  transcription  factor  and  acts 
as  a  transcription  repressor.  The  effect  of  NFKBIE  upre¬ 
gulation  on  apoptosis  is  unclear  because  it  was  reported 
to  exhibit  a  cell  line  dependence  [37]. 

Cell  division  cycle  25  homolog  A,  B,  and  C  expression 
showed  a  remarkable  regulation;  e.g.,  Cdc25,  which 
plays  a  key  role  in  promoting  progress  through  S  phase 
[38],  showed  1.7  FC,  while  Cdc  25  B  and  C,  which  regu¬ 
late  and  promote  entry  into  mitosis  [39,40],  exhibited 
-4.3  and  -1.92  FC,  respectively.  Similarly,  cyclin  -depen¬ 
dent  kinase  regulators,  cyclin  E2  and  cyclin  Dl,  which 
are  strongly  tied  to  a  variety  of  cancers  [41-43],  showed 
upregulation  of  3.2  and  2.2  FC,  respectively. 

Other  cell  cycle  regulators  that  have  overlapping  func¬ 
tions  with  apoptosis  mechanisms  showed  sharp  downre- 
gulation.  This  include  members  of  the  p21  RAS  group 
such  as  KRAS  (v-I<i-ras2  Kirsten  rat  sarcoma  viral  onco¬ 
gene  homolog)  and  Ras  homolog  or  the  RHO  GTPase 
protein  group  including  RHOB,  RHOQ,  RHOT1,  and 
RND3  (Fig.  6). 

Ephrine  signaling  pathway  (ESP) 

The  largest  group  of  receptor  tyrosine  kinases  (RTKs), 
Eph  receptor  family  and  their  ligands,  Ephrines,  were 
identified  among  the  main  pathways  that  were  influ¬ 
enced  by  MPV  infection.  Genes  in  the  7  hpi  data  set 
were  enriched  to  ESP  at  ratio  of  13.8%  of  the  total  195 
pathway  genes  with  a  P-value  of  1.37e-4  (Fig.  5).  The 
percentage  of  pathway  genes  regulated  by  infection  was 
lower  during  the  3  hpi  time  point,  reaching  only  a  ratio 
of  10.3%  with  P-value  of  5.95e-5,  and  ranking  fourth 
among  influenced  pathways  by  MPV  infection.  Trans¬ 
membrane  receptors  Eph  A  class,  which  contains  eight 
members  of  ephrine  receptors  and  represent  one  of  two 
classes  forming  the  entire  receptors  family  [44],  exhib¬ 
ited  significant  suppression  associated  with  downregula- 
tion  of  key  kinases  and  protein  groups  implicated  in  cell 
morphology,  cell  proliferation,  and  cell  repulsion.  Of 
these  molecules,  we  mention  FYN  oncogene,  RAS  p21 
protein  activator  (GTPase  activating  protein  RAS-GAP) 
1,  p21  protein  (Cdc42/Rac) -activated  kinase  (PAI<),  and 
Rho-associated,  coiled-coil  containing  protein  kinase  1 
(member  of  ROCK  group)  (Fig.  7).  Interestingly,  the 
other  class  of  ephrines  family,  Eph  B,  didn’t  show  any 
significant  change  in  its  expression  in  both  time  points. 
However  many  molecules  downstream  of  Eph  B  in  the 
signal  cascade  were  downregulated,  including  intersectin 
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Figure  5  Major  influenced  canonical  pathways  of  MPV-infected  MK2  cells.  Relevant  pathways  are  shown  at  3  (A)  and  7  (B)  hpi.  The  primary 
Y-axis  shows  the  -log  (P-value)  of  the  probability  for  genes  in  data  set  to  associate  with  identified  pathway  by  chance.  A  threshold  P-value  of 
0.05  or  a  -log  (P)  =  1.3  is  presented  in  dotted  line.  The  ratio  of  the  number  of  genes  from  the  data  set  that  map  to  given  pathway  divided  by 
the  total  number  of  genes  that  map  to  the  canonical  pathway  is  shown  in  solid  line.  All  enlisted  pathways  met  ratios  >  0.1  in  at  least  one  of  the 
two  time  points  (secondary  Y-axis). 

..  j 
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Figure  6  Heat-map  of  differentially  expressed  genes  in  MPV 
infected  MK2  cells  that  mapped  to  molecular  mechanisms  of 
cancer  pathway.  Columns  represent  expression  at  3  and  7  hpi  in  a 
triplicate  (A,  B,  C).  Each  row  represents  one  gene  that  met  analysis 
cutoff  of  average  fold  change  >1.8  in  at  least  one  of  the  two  time 
points  and  P-value  <  0.05  in  both.  A  gradient  of  green  and  red 
colors  represent  low  and  high  relative  fold  change  of  gene 
expression  to  mock-  infected  cells.  Chance  for  random  association 
of  listed  genes  at  3  or  7  hpi  in  this  pathway  is  2.82e-6  or  1.25e-5, 
respectively. 


1  (SH3  domain  protein),  c-abl  oncogene  1-RTK  (ABL), 
and  Wiskott-Aldrich  syndrome-like  (WASL),  cofilin  2 
(CFL2),  actin-related  protein  2  homolog  (ARP2),  and 
mitogen-activated  protein  kinase  4  (NcK)  (Fig.  7). 

Two  proteins  were  up-regulated  in  this  pathway, 
ARPC5L  (actin-related  protein  2/3  complex,  subunit  5- 
like),  which  functions  as  an  actin-binding  protein  and  is 
involved  in  regulation  of  actin  filament  polymerization 
[45]  and  GNAS  complex  locus,  which  is  a  member  of 
G -protein  alpha  subunit  that  has  an  indirect  effect  on 
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Figure  7  Cluster  Analysis  of  expression  profile  for  genes  of 
ephrin  receptor  signaling  pathway  influenced  by  MPV 
infection.  Columns  of  the  heat-map  represent  expression  of  genes 
transcripts  at  3  and  7  hpi  in  a  triplicate  (A,  B,  C).  Each  row 
represents  one  gene  that  met  the  analysis  parameters  of  si. 8 
average  fold  change  cutoff  in  at  least  one  of  the  two  time  points 
and  P-value  <  0.05  in  both.  Gradient  green  and  red  color  represents 
low  and  high  relative  expression  fold  change  to  mock  infected  cells, 
respectively.  Chance  for  random  association  of  listed  genes  at  3  or  7 

hpi  with  this  pathway  is  5.95e-5  or  1 ,37e-4,  respectively, 
v _ _ _ _ _ 


signaling  of  ephrine  pathway.  The  importance  of  this 
pathway  is  due  to  its  involvement  in  the  modulation  of 
integrin  activity,  actin  reorganization,  and  cellular 
dynamics,  which  are  all  central  to  viral  motility,  exit  and 
invasion  of  proximal  cells. 

Ataxia  telangiectasia  mutated  protein  signaling  pathway 

Ataxia  telangiectasia  mutated  protein  (ATM)  is  a  cell 
kinase  that  phosphorylates  a  wide  range  of  substrates  in 
many  pathways  regulating  cell  cycle,  DNA  repair,  apop¬ 
tosis,  and  cell  survival.  Data  sets  analysis  showed  mild 
ATM  downregulation  of  1.3  and  1.5  FC  in  time  points  3 
and  7  hpi,  respectively.  Nonetheless,  ATM  pathway  was 
recognized  in  our  analysis  as  one  of  the  major  pathways 
influenced  by  MPV  infection.  Out  of  the  52  genes  com¬ 
prised  in  the  ATM  pathway,  eight  genes  exhibited  sig¬ 
nificant  expression  modulation  with  P-value  of  1.87e-3 
at  3  hpi.  The  number  of  genes  regulated  at  7  hpi 
increased  to  10  genes,  resembling  19.2%  of  the  total 
pathway  components  with  P-value  of  4.53e-3.  The 
observed  suppression  of  ATM  results  in  weaker 
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signaling  of  the  pathway  and  consequently  reduces 
response  to  DNA  damage  and  routine  DNA  repair, 
additionally,  it  disables  a  major  apoptosis  activation  pro¬ 
cess.  The  MPV  infection  influence  on  the  pathway  is 
exacerbated  further  by  the  downregulation  of  many  vital 
ATM  substrates  needed  for  proper  ATM  signaling.  This 
includes  activating  transcription  factor  2  (ATF2),  c-abl 
oncogene  1  receptor  tyrosine  kinase  (c-Abl),  and  breast 
cancer  1  (BRCA1)  (Fig.  8). 

Pathways  involved  in  cell  cycle  regulation 
A  significant  proportion  of  the  genes  in  data  sets  3  and 
7  hpi  were  enriched  in  three  major  cell  cycle  regulation 
pathways.  Out  of  the  36  molecules  found  in  cell-cycle 
regulation  by  BTG  family  pathway,  seven  molecules 
were  found  in  each  of  the  3  and  7  hpi  data  sets  with 
P-values  of  8.15e-4  at  3  hpi  and  1.48e-2  at  7  hpi  indicat¬ 
ing  the  involvement  of  this  pathway  in  early  time  point. 
Essential  genes  to  pathway  signaling  and  G1  arrest  were 
severely  suppressed  (Figure.  9A).  Retinoblastoma  1  (Rb) 
and  B-cell  translocation  gene-1  (BTG1)  were  downregu- 
lated  in  both  time  points  and  reached  -9.26  and  -15.86 
FC  by  7  hpi.  Furthermore,  cyclin  D1  and  cyclin  E  exhib¬ 
ited  4.18-  and  2.77-fold  increase  at  3  hpi  respectively, 
and  remained  upregulated  to  similar  levels  at  7  hpi  time 
point.  The  products  of  these  two  genes  function  as  reg¬ 
ulators  that  inactivate  the  growth  suppression  activity  of 
Rb  indirectly  by  promoting  CDI<  kinases  [46-48]. 


Phosphorylation  of  Rb  by  CDI<  kinase  inactivates  its 
function,  leading  to  the  release  of  E2F  and  cell  cycle 
progression  [49,50].  The  net  outcome  of  this  pattern  of 
gene  expression  in  this  pathway  is  cell  cycle  release 
from  potential  arrest  in  G1  phase  and  block  of  BTG1- 
induced  apoptosis. 

The  second  pathway  contributing  to  cell  cycle  regula¬ 
tion  that  was  influenced  by  infection  is  the  cell  cycle 
Gl/S  checkpoint.  The  pathway  consists  of  59  molecules 
from  which  eight  molecules  or  13.6%  showed  significant 
regulation  in  both  data  sets  (Fig.  9B).  The  calculated  P- 
values  for  random  implication  of  this  pathway  were 
2.71e-3  and  5.02e-2  for  the  3  and  7  hpi,  respectively. 
Similar  to  regulation  by  the  BTG  protein  family  path¬ 
way,  phosphorylation  of  retinoblastoma  Rb  by  CDI<2 
and  CDK4/6  induces  the  release  of  transcription  factor 
E2F  from  an  inhibitory  complex,  which  in  turn  pro¬ 
motes  the  transcription  of  necessary  molecules  for  Gl/S 
phase  progression  [50-52],  Cell-division-cycle  25  homo¬ 
log  A  (Cdc25A),  which  activates  CDI<2  and  CDK4/6 
[53-55],  showed  strong  upregulation  suggesting 
enhancement  of  kinase  activity  that  promotes  Gl/S 
phase  progression.  Transforming  growth  factor  |3 
(TGFp)  that  stimulates  Gl/S  check  point  exhibited 
almost  steady  suppression  of  -1.48  and  -1.52  FC  in  both 
time  points  [56,57].  Another  mechanism  of  Gl/S  check¬ 
point  control  is  degradation  of  cyclin  D  and  E  after 
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Figure  8  Cluster  analysis  of  expression  profile  for  genes  of  ataxia  telangiectasia  mutated  (ATM)  pathway  Columns  of  the  heat-map 
represent  expression  of  genes  transcripts  at  3  and  7  hpi  time  points  in  a  triplicate  (A,  B,  C).  Each  row  represents  one  gene  that  met  the  analysis 
criteria  of  >1.8  average  fold  change  cutoff  in  at  least  one  of  the  two  time  points  and  P-value  <  0.05  in  both.  Gradient  green  and  red  color 
represents  low  and  high  relative  expression  fold  change  to  mock  infected  cells,  respectively.  Chance  for  random  association  of  listed  genes  at  3 
or  7  hpi  with  this  pathway  is  1 .87e-3  or  4.53e-3,  respectively. 


Alkhalil  et  al.  Virology  Journal  2010,  7:173 
http://www.virologyj.eom/content/7/1/173 


Page  14  of  19 


oo  on  on  t'--  t-' 

<  co  kJ  <  go  GJ 


A. 


ABL1 

RBI 

E2F5 

RBL2 

CDKN2B 

CDC25A 

CCNDl^ 

CCNE2 


> 


Gl/S  Checkpoint 
Regulation 


B. 


CDK7 

GADD45A 

CDC25C 

RPS6KA1 

CDC25B 

PRKDC 

CCNB2 

T0P2A 

SKP1 

T0P2B 

BRCA1 


G2/M  Checkpoint 
Regulation 


C. 


BTG1 

PPP2R5E 

RBI 

PPP2R5C 

E2F5 

CCNE2 

CCND1 


Regulation  by  BTG 


Family  Proteins  Pathway 


<  -2  0  +2  < 


down  regulated  up  regulated 


Figure  9  Cluster  analysis  of  differentially-  expressed  genes  involved  in  cell  cycle  regulation  Heat-map  of  gene  expression  levels  at  3  and 
7  hpi  in  triplicate  (A,  B,  C).  Folds  of  change  in  gene  expression  represented  as  a  gradient  of  green  and  red  color  for  low-  and  high-expression 
intensity,  respectively.  Each  row  represents  one  gene  that  met  analysis  cutoff  of  average  fold  change  >1.8  in  at  least  one  of  the  two  time  points 
and  P-value  s  0.05  in  both.  Genes  fell  in  three  cell  cycle  pathways  with  some  genes  functioning  in  more  than  one  pathway  (A)  cell  cycle:  Gl/S 
checkpoint  regulation  pathway  (B)  cell  cycle:  G2/M  checkpoint  regulation  pathway  (C)  cell  cycle  regulation  by  BTG  family  proteins  pathway. 


activation  by  glycogen  synthase  kinase  |3  (GSI<P)  [58,59]. 
Both  time  points  show  upregulation  of  cyclin  D  and  E 
expression  with  no  change  in  GSI<P  expression  despite 
mild  reduction  in  TGFp.  This  multigene  expression  pat¬ 
tern  in  both  data  sets  points  to  destabilized  Gl/M  con¬ 
trol  checkpoint  and  enhancement  of  cell  cycle 
progression  through  this  phase,  especially  in  the  early 
time  point  as  suggested  by  lower  P-value  and  greater 
gene  regulation. 

The  third  cell-cycle-related  pathway  identified  in  our 
analysis  was  the  G2/M  DNA  damage  checkpoint 


regulation  pathway.  About  14%  of  the  total  43  molecules 
comprised  in  this  pathway  exhibited  significant  expres¬ 
sion  modulation  upon  infection  at  time  point  3  hpi  with 
P-value  of  5.02e-2.  Time  point  7  hpi  showed  regulation 
of  20.9%  of  total  pathway  with  P-value  of  1.19e-3,  imply¬ 
ing  a  wider  pathway  role  during  later  stage  of  infection. 
Progression  of  cell  cycle  into  mitosis  requires  the  activa¬ 
tion  of  cell-division  cycle-2  (Cdc2)  by  the  rapid  depho¬ 
sphorylation  of  its  tyrosine  15  (Tyr-15)  via  action  of 
Cdc25  and  subsequent  binding  to  cyclin  B  [60,61]. 
Damaged  DNA  during  G2  phase  triggers  ATM  and 
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ataxia  telangiectasia  mutated  and  rad3-related  (ATR) 
[62-64],  which  in  turn  activates  two  key  kinases,  Chk2 
and  Chkl [65,66].  The  latter  kinase  inhibits  Cdc25B  and 
Cdc25C  which  function  as  a  phosphorylase  that  activate 
Cdc2  before  binding  to  cyclin  B  and  entry  into  mitosis 
[67].  Our  results  showed  that  ATR/ATM,  and  Chkl, 
which  activate  sequentially  in  response  to  DNA  abnorm¬ 
alities,  exhibited  no  or  only  a  subtle  change  in  their 
expression  in  both  time  points,  suggesting  retention  of 
their  Cdc25B/C  inhibitory  effect  as  in  non-infected  cells. 
Three  essential  molecules  needed  for  entry  into  mitosis 
exhibited  substantial  downregulation;  Cdc25B/C  and 
cyclin  B  exhibited  -4.33,  -1.91,  -4.48  FC  at  7  hpi  and 
-1.81,  -2.27,  -1.71  FC  at  3  hpi  respectively.  The  observed 
regulation  of  these  genes  would  intensify  this  checkpoint 
signaling  leading  to  a  further  delay  in  entry  into  mitosis 
(Fig.  9C). 

Discussion 

Apoptosis  is  a  natural  controlled  cell  death  mechanism 
triggered  by  diverse  stimuli  to  maintain  tissue  homeosta¬ 
sis  and  eliminate  abnormal  or  infected  cells  [68].  Because 
apoptosis  represents  an  important  part  of  antiviral  host 
response,  poxviruses  developed  numerous  ways  to  target 
it  and  disrupt  its  function  [69].  Diverse  anti-apoptotic 
viral  strategies  are  identified  in  different  poxviruses,  e.g., 
Molluscum  contagiosum  virus  (MCV)  inhibits  caspase-8 
by  expressing  MC159  gene,  which  encodes  a  protein  that 
binds  procaspase-8  and  Fas-associated  death  domain, 
thereby  inhibiting  death  receptor-induced  apoptosis 
mediated  by  Fas,  TNF,  or  TRAIL  receptors  [70].  MC066 
gene  encodes  a  protein  with  glutathione  peroxidase-like 
function  to  convert  oxygen-reactive  species  to  neutral 
molecules,  hence  preventing  apoptosis  triggered  by 
increased  oxidative  stress  associated  with  infection  [71]. 
Little  is  known  about  apoptosis  in  cells  infected  with 
MPV,  but  related  orthopoxviruses  exhibit  clear  anti- 
apoptotic  functions.  Cowpox  virus  for  instance  expresses 
a  protein  that  can  block  apoptosis  in  multiple  ways.  One 
of  the  most  potent  anti-apoptotic  proteins  is  the  cytokine 
response  modifier  A  (CrmA)  which  inhibits  caspase-8, 
caspase-10,  and  blocks  garnzyme  B-mediated  apoptosis 
[72,73].  Similarly,  Vaccinia  viruses  use  SPI-2  family  pro¬ 
tein  member  of  serine  protease  inhibitors  (serpins) 
encoded  by  B13R  gene  to  block  apoptosis  induced  by 
death  receptors  [74].  The  same  virus  inhibits  apoptosis 
induced  by  RNA-dependent  protein  kinase  (PI<R)  using 
specific  PI<R  inhibitors  encoded  by  E3L  and  I<3L  genes 

[75] ,  and  apoptosis  induced  by  loss  of  mitochondria 
outer  membrane  potential  by  expression  of  F1L  gene 

[76] , 

The  presence  of  many  viral  proteins  that  block  apop¬ 
tosis  at  multiple  points  suggests  that  apoptosis  is  detri¬ 
mental  to  viral  survival.  However,  the  observed 


downregulation  of  Bcl-2,  PUMA,  and  PAI<2  and  upregu- 
lation  of  NOXA  and  caspase-3  with  the  negligible 
change  in  other  apoptotic  genes  in  our  data  is  more 
consistent  with  apoptosis  induction.  Although  many 
pathogens  and  viruses  other  than  poxviruses  are 
reported  to  promote  apoptosis  [77],  it  is  unlikely  that 
MPV  will  differ  from  other  poxviruses  in  their  common 
anti-apoptotic  trend,  and  the  observed  divergence 
between  the  regulation  of  apoptosis-specific  genes  in 
MPV-infected  cells  and  overall  anti-apoptotic  outcome 
seen  in  other  poxvirus-infected  cells  suggests  an  anti- 
apoptotic  viral  mechanism  that  functions  downstream  of 
apoptosis  induction  in  the  host.  MPV  genes  involved  in 
blocking  apoptosis  remain  unknown,  but  our  data  sug¬ 
gests  the  presence  of  an  ortholog(s)  of  Vaccinia  virus 
(F1L)  gene  in  MPV,  which  acts  directly  on  the  mito¬ 
chondria  in  Bcl-2  like  manner. 

Double-stranded  DNA  mammalian  viruses  have  large 
genomes  reaching  289  Kbp  or  few  hundreds  of  micro¬ 
meters  in  length  as  in  the  case  of  fowlpox  virus  [78]. 
Compacting  the  genome  is  an  indispensible  biological 
task  due  to  the  rich  negative  charge  and  large  size  of 
DNA  molecules.  In  eukaryotes,  this  was  solved  by  wrap¬ 
ping  DNA  around  a  heterooctameric  molecule  com¬ 
posed  of  four  dimerized,  positively  charged  core 
histones  to  form  nucleosomes.  Further  stability  is 
brought  about  by  linker  histones  that  bind  DNA  later¬ 
ally.  Nucleosomes,  in  a  nucleofilaments  form  or  in  a 
form  of  higher-order-structure  assemblies,  chromatin, 
play  an  essential  role  in  the  regulation  of  gene  expres¬ 
sion  and  chromatin  remodeling  via  acetylation  and  dea¬ 
cetylation  of  the  N-terminal  histone  tails  protruding 
from  nucleosomes.  Viruses  exhibit  similar  architecture 
in  their  DNA  compaction.  Staining  vaccinia-infected 
cells  with  osmium  ammine-S02  revealed  areas  with  con¬ 
centrated  viral  DNA  in  two  predominant  DNA  config¬ 
urations  varying  depending  on  the  level  of  viral 
synthetic  activities.  Encapsidated  genomes  exhibited  a 
nucleosome  structure  like  those  observed  in  resting 
eukaryotes,  which  was  in  agreement  with  biochemical 
data  showing  the  supercoiled  organization  of  nucleopro- 
teins  extracted  from  Vaccinia  viruses  [79].  On  the  other 
hand,  like  active  cellular  chromatin,  non-encapsidated 
viral  genomes  exhibited  extended  DNA  features  [80]. 
Similar  variation  in  DNA  density  was  observed  in 
Herpes  simplex  viruses  (HSV),  another  double  stranded 
DNA  virus,  with  two  well  described  configurations  of 
euchromatin  or  heterochromatin  that  correlated  strongly 
with  lytic  or  latent  HSV  infection  stages,  respectively. 
Generally,  little  is  known  about  chromatin  potential  in 
poxviruses  as  an  antiviral  mechanism.  Our  results 
showed  an  interesting  downregulation  of  five  essential 
histone  expression  regulation  factors  and  five  enzymes 
regulating  chromatin  dynamics  in  MPV-infected  cells.  It 
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is  unclear  if  these  results  are  part  of  the  host  response, 
which  include  chromatin-mediated  silencing  of  the  viral 
genome  and  activation  of  DNA  damage  [81],  or  part  of 
the  viral  strategies  to  take  over  its  host.  However,  our 
observation  predicts  an  important  role  for  histone 
expression,  histone  posttranslational  modification,  and 
dynamic  exchanges  of  chromatin  in  host-poxvirus  inter¬ 
actions.  Recent  work  suggested  a  role  for  the  viral  A32L 
gene  of  Vaccinia  virus  in  DNA  packaging  based  on 
sequence  similarities  with  the  product  of  gene  I  of  fila¬ 
mentous  single-stranded  DNA  bacteriophages  and  the 
Iva2  gene  of  adenoviruses.  Both  of  these  genes  are 
ATPases  involved  in  DNA  packaging.  Additional  Vacci¬ 
nia  genes  that  map  to  16  or  II  telomere-binding  pro¬ 
teins  are  believed  to  play  roles  in  DNA  packaging, 
because  mutants  of  VACV  in  either  of  these  two  genes 
fail  to  exhibit  normal  DNA  packaging  at  different  mor¬ 
phogenesis  stages.  The  sharp  upregulation  we  observed 
in  three  out  the  four  core  histones,  and  the  striking 
similarities  in  DNA  compaction  architecture  observed  in 
eukaryotes  and  some  viruses,  with  absence  of  known 
poxvirus  proteins  that  exhibit  histone-like  properties 
makes  it  tempting  to  hypothesize  a  role  of  host  cell 
histones  in  viral  DNA  compaction  and  nucleosome 
formation,  especially  that  similar  involvement  was 
described  recently  in  simian  virus  40  (SV40)  DNA  com¬ 
paction  [82] 

Our  analysis  identified  ephrin  receptor  pathway  (ERP) 
as  a  major  influenced  pathway  in  infected  cells.  This 
might  be  due  to  either  increased  cell  to  cell  communica¬ 
tion  by  signaling  through  this  receptor  tyrosine  kinases 
(RTI<)  family  in  response  to  infection,  or  to  the  presence 
of  many  pleiotropic  genes  that  are  found  in  ERP,  and 
simultaneously  have  essential  roles  in  cytoskeleton  reor¬ 
ganization  or  actin  polymerization.  Intracellular  viral 
motility  and  morphogenesis  of  Vaccinia  virus  into  a 
cell-associated  enveloped  virion  (CEV)  and  extracellular 
enveloped  virions  (EEV)  forms  are  shown  to  be  driven 
by  interactions  of  host  microtubules  and  Vaccinia  A27L, 
A17L  and  A14L  genes.  Furthermore,  egress  of  Vaccinia 
particles  in  EEV  form  and  direct  cell-to-cell  virus  disse¬ 
mination  is  propelled  by  actin  tail  formation,  which 
involves  the  interactions  of  Vaccinia  transport  genes, 
including  A36R,  F12L,  and  host  proteins  such  as  Src 
family  kinases  (SFI<),  Nek,  WIP,  N-WASP,  Arp  2/3  to 
promote  actin  filaments  nucleation.  Actin  polymeriza¬ 
tion  produces  microvilli  at  cell  surface  that  lift  CEV  and 
project  it  on  adjacent  cells  to  finally  deliver  the  virus 
with  minimal  exposure  to  host  immune  system.  Our 
results  confirm  regulation  of  many  principal  signaling 
components  involved  in  actin  cytoskeletal  dynamics,  and 
introduce  additional  infection  regulated  genes  with  func¬ 
tions  related  to  microtubules  signaling.  This  includes 
intersectin  1  (SH3  domain  protein)  gene,  which  encodes 


a  cytoplasmic  membrane-associated  protein  that  indir¬ 
ectly  coordinates  endocytic  membrane  traffic  with  the 
actin  assembly  machinery,  Rho-effector  ROCK1  serve  a 
number  of  key  cellular  functions,  such  as  morphological 
differentiation  and  cell  motility  which  are  closely  asso¬ 
ciated  with  changes  in  cytoskeletal  dynamics  [83].  Addi¬ 
tionally,  RAS  p21  protein  activator  (GTPase  activating 
protein)  [84],  v-I<i-ras2  Kirsten  rat  sarcoma  viral  onco¬ 
gene  homolog  [85]  and  SOS2  [86]  are  crucial  genes  in 
polymerization  of  actin  filaments  and  cytoskeleton 
reorganization. 

Ion  channels  represent  an  intriguing  and  novel  class  of 
genes  that  were  impacted  by  MPV  infection.  We  identi¬ 
fied  10  genes  encoding  nine  ion  channels  and  a  trans¬ 
porter  that  underwent  increasing  suppression  during 
infection.  Most  of  these  channels  localize  to  cell  mem¬ 
brane,  and  collectively  contribute  to  transport  of  all 
essential  ions  involved  in  maintenance  of  cell  membrane 
potential  and  osmolarity  homeostasis.  While  mechan¬ 
isms  of  transport  modulation  have  been  described  pre¬ 
viously,  as  in  the  indirect  consequences  of  Ras,  Rho,  and 
Rab  small  GTPases  regulation  [87],  its  effect  on  viral 
infections  and  global  cell  biology  remain  unclear  except 
for  a  recent  report  describing  the  interaction  of  myxoma 
poxvirus  protein  M11L  with  mitochondrial  permeability 
transition  pore  and  its  role  in  delaying  apoptosis  in  host 
cells  [88].  The  downregulation  trend  of  channel  expres¬ 
sion  identified  here  pose  many  intriguing  questions, 
especially  in  the  light  of  evolving  evidence  in  support  of 
ion  channels  role  in  virus  release  [89]  and  infected  cells 
rupture  [90]. 

Progression  of  the  cell  cycle  is  tightly  regulated  pro¬ 
cess  with  many  redundant  checkpoints  that  ensure 
proper  transition  across  cell  cycle  phases.  Our  results 
showed  significant  modulation  in  the  expression  of 
many  genes  that  play  essential  roles  in  cell  cycle  regula¬ 
tion,  which  led  to  the  identification  of  ATM  signaling, 
G2/M  DNA  damage  checkpoint,  regulation  by  BTG 
family  protein,  and  Gl/S  checkpoint  as  major  influenced 
pathways  during  MPV  infection.  A  core  cell  cycle  regu¬ 
lation  gene,  Cdc25  kinase,  has  three  essential  homologs 
Cdc  25A/B/C  that  exhibited  significant  regulation  upon 
MPV  infection.  While  Cdc25B/C  showed  downregula¬ 
tion  favoring  cell  arrest  in  G2  phase,  Cdc2A  exhibited 
upregulation,  favoring  S  phase  progression.  The  impact 
of  this  mode  of  cell  cycle  regulation  on  viral  infection 
remains  unknown.  However,  cell  arrest  in  G2  phase  was 
described  in  other  viral  infections  including  human 
immunodeficiency  virus  (HIV)  and  was  found  to  be 
mainly  mediated  by  viral  protein  R  (Vpr)  [91-93].  While 
many  of  the  G2/M  DNA  damage  checkpoint  pathway 
genes  are  known  to  be  modulated  during  HIV  infection, 
Vpr  seems  to  induce  cell  arrest  by  molecular  mechan¬ 
isms  other  than  the  classic  DNA  check  point  [94]. 
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Recently,  evidence  supporting  a  role  for  PP2A  in  Vpr- 
induced  arrest  has  emerged,  and  was  substantiated 
further  by  other  studies  in  support  of  PP2A  being  a 
common  target  during  infection  with  other  viruses, 
including  simian  virus  40  (SV40),  polyoma  virus,  human 
T  lymphotrophic  retrovirus  and  adenovirus  [95].  Our 
results  showed  significant  downregulation  in  two  PP2A 
isoforms,  regulatory  subunit  B’  gamma  isoform 
(PPP2R5C)  and  protein  phosphatase  2,  regulatory  subu¬ 
nit  B’  epsilon  isoform  (PPP2R5E),  suggesting  that  the 
induction  mechanisms  of  G2  arrest  in  MPV  infection 
might  be  similar  to  those  observed  in  other  viruses. 
Because  genetically  diverse  viruses  seem  to  induce  the 
same  G2  arrest  response  in  different  infected  cells,  it  is 
likely  that  this  response  has  an  important  function  and 
might  be  part  of  antiviral  host  defenses.  While  some  of 
the  viral  genes  eliciting  this  response  are  being  identified 
as  Vpr  in  HIV,  and  E4  of  F4  and  HTLV  tax  protein  in 
adenoviruses,  the  MPV  gene  inducing  this  response 
remains  unknown.  Other  important  genes  in  cell  cycle 
regulation  showing  expression  favoring  progression  of 
cell  cycle  and  arrest  only  in  G2  phase  include  Rb,  E2F, 
cyclins,  BTG1,  and  BRCA1. 

In  this  study  we  combined  microarray  with  data 
mining  and  statistical  analysis  to  identify  important 
interfaces  of  host-pathogen  interaction.  Our  results 
aligned  nicely  with  previous  reports  carried  out  using 
viruses  from  the  same  or  different  genus,  and  provided 
new  set  of  genes  that  play  important  roles  in  MPV 
infection.  Further  work  is  warranted  to  validate  and 
examine  the  potential  of  these  genes  in  antiviral 
therapies. 

Summary 

Using  microarrays,  we  studied  MPV-induced  changes  in 
gene  expression  of  Macaca  mulatta  kidney  epithelial 
cells  to  identify  major  host-virus  interaction  interfaces. 
Infection  stimulated  a  marked  modulation  in  some 
2,702  host  genes  representing  ~  5.7%  of  total  interro¬ 
gated  host  transcripts.  The  majority  of  genes  (89.08%) 
underwent  1.5-folds  downregulation  or  more.  While 
downregulated  genes  exhibited  a  steady  trend  during  the 
study,  upregulated  genes  showed  more  time-dependent 
regulation  intensity.  Further  data  analysis  showed  that 
regulated  genes  cluster  into  distinctive  functional 
classes,  canonical  pathways  and  networks  that  can  be 
linked  to  established  viral  biogenesis.  Our  results  intro¬ 
duce  a  set  of  host  genes  and  novel  pathways  for  further 
evaluation  as  targets  for  potential  use  in  developing  new 
antiviral  therapies. 

Abbreviations 

CITED2:  Cbp/p300-i interacting  transactivator,  with  Glu/Asp-rich  carboxy- 
terminal  domain,  2;  NC0A3:  nuclear  receptor  coactivator  3;  CREB:  cAMP 


responsive  element  binding  protein  1;YY1:  YY1  transcription  factor;  HDAC2: 
histone  deacetylase  2;  FBX011:  F-box  protein  11;  PRMT3:  protein  arginine 
methyltransferase  3;  MYST2:  MYST  histone  acetyltransferase  2;  MYCBP2:  MYC 
binding  protein  2;  RARS2:  arginyl-tRNA  synthetase  2,  mitochondrial;  NOXA: 
NADPH  oxidase  activator;  Bcl2:  B-cell  CLL/lymphoma  2;  BBC3  or  PUMA:  BCL2 
binding  component  3;  PAK2:  p21  protein  (Cdc42/Rac)-activated  kinase  2; 
Cdc25:  cell  division  cycle  25;  KRAS:  v-Ki-ras2  Kirsten  rat  sarcoma  viral 
oncogene  homolog;  RHOB:  ras  homolog  gene  family,  member  B;  RHOQ:  ras 
homolog  gene  family,  member  Q;  RH0T1:  ras  homolog  gene  family, 
member  T1;  RND3:  Rho  family  GTPase  3;  FYN:  FYN  oncogene  related  to  SRC, 
FGR,  YES;  ROCK:  Rho  Associated  Kinase,  Rho  kinase;  ABL1/2:  c-abl  oncogene 
1  or  2,  receptor  tyrosine  kinase;  WSAL:  Wiskott-Aldrich  syndrome-like;  CFL2: 
cofilin  2;  ARP2/3:  actin  related  protein  2/3;  ARPC5L:  actin-related  protein  2/3 
complex,  subunit  5-like;  GNAS:  GNAS  complex  locus;  ATF2:  activating 
transcription  factor  2;  c-Abl:  c-abl  oncogene  1  receptor  tyrosine  kinase; 
BRCA1:  breast  cancer  1;  Rbl:  Retinoblastoma  1;  BTG1:  B-cell  translocation 
gene-1;  CDK:  cyclin-dependent  kinase;  TGF0:  Transforming  growth  factor  p; 
GSKp:  glycogen  synthase  kinase  p;  ATR:  ataxia  telangiectasia  mutated  and 
rad3-related;  ChK:  choline  kinase;  TRAIL  or  TNFSF10:  tumor  necrosis  factor 
(ligand)  superfamily,  member  10;  FAS:  Fas  (TNF  receptor  superfamily, 
member  6);  CrmA:  cytokine  response  modifier  A;  SFK:  Src  family  kinases; 
S0S2;  son  of  sevenless  homolog  2;  PP2A:  Protein  Phosphatase  Type2a;  HTLV: 
human  T-cell  leukemia  virus 

Acknowledgements 

Opinions,  interpretations,  conclusions,  and  recommendations  are  those  of 
the  author  and  are  not  necessarily  endorsed  by  the  U.S.  Army. 

Microarray  data  have  been  submitted  to  the  Gene  Expression  Omnibus 
(GEO)  and  can  be  searched  using  the  record  ID:  GSE21001. 

Author  details 

United  States  Army  Medical  Research  Institute  of  Infectious  Diseases,  Fort 
Detrick,  Maryland  21702,  USA.  2Division  of  Pathology,  Walter  Reed  Army 
Institute  of  Research,  Silver  Spring,  MD,  20910,  USA. 

Authors'  contributions 

AA  was  responsible  for  design,  conduct,  and  completion  of  this  work,  as 
well  as  for  data  analysis  and  writing  of  this  manuscript.  RH  and  MJ  were 
instrumental  in  data  processing  and  statistical  analysis  of  microarray  data.  JH 
and  MA  contributed  to  microarray  data  validation  using  RT  PCR.  SI  was  the 
Principal  Investigator  and  is  primarily  responsible  for  all  aspects  of  the 
funding,  research  design,  interpretation,  and  writing  of  this  manuscript.  All 
authors  read  and  approved  the  final  manuscript. 

Competing  interests 

The  authors  declare  that  they  have  no  competing  interests. 

Received:  7  May  2010  Accepted:  28  July  2010  Published:  28  July  2010 

References 

1.  Breman  JG,  Kalisa  R,  Steniowski  MV,  Zanotto  E,  Gromyko  Al,  Arita  I:  Human 
monkeypox,  1970-79.  Bull  World  Health  Organ  1980,  58:165-182. 

2.  Likos  AM,  Sammons  SA,  Olson  VA,  Frace  AM,  Li  Y,  Olsen-Rasmussen  M, 
Davidson  W,  Galloway  R,  Khristova  ML,  Reynolds  MG,  ef  ah  A  tale  of  two 
clades:  monkeypox  viruses.  J  Gen  Virol  2005,  86:2661-2672. 

3.  Chen  N,  Li  G,  Liszewski  MK,  Atkinson  JP,  Jahrling  PB,  Feng  Z,  Schriewer  J, 
Buck  C,  Wang  C,  Lefkowitz  EJ,  ef  al:  Virulence  differences  between 
monkeypox  virus  isolates  from  West  Africa  and  the  Congo  basin. 

Virology  2005,  340:46-63. 

4.  Shchelkunov  SN,  Totmenin  AV,  Safronov  PF,  Mikheev  MV,  Gutorov  W, 
Ryazankina  01,  Petrov  NA,  Babkin  IV,  Uvarova  EA,  Sandakhchiev  LS,  ef  al: 

Analysis  of  the  monkeypox  virus  genome.  Virology  2002,  297:172-194. 

5.  Carter  GC,  Law  M,  Hollinshead  M,  Smith  GL:  Entry  of  the  vaccinia  virus 
intracellular  mature  virion  and  its  interactions  with  glycosaminoglycans. 

J  Gen  Virol  2005,  86:1279-1290. 

6.  Hsiao  JC,  Chung  CS,  Chang  W:  Vaccinia  virus  envelope  D8L  protein  binds 
to  cell  surface  chondroitin  sulfate  and  mediates  the  adsorption  of 
intracellular  mature  virions  to  cells.  J  Virol  1999,  73:8750-8761. 

7.  Chung  CS,  Hsiao  JC,  Chang  YS,  Chang  W:  A27L  protein  mediates  vaccinia 
virus  interaction  with  cell  surface  heparan  sulfate.  J  Virol  1998, 
72:1577-1585. 


Alkhalil  et  al.  Virology  Journal  2010,  7:173 
http://www.virologyj.eom/content/7/1  /1 73 


Page  18  of  19 


8.  Lin  CL,  Chung  CS,  Heine  HG,  Chang  W:  Vaccinia  virus  envelope  H3L 
protein  binds  to  cell  surface  heparan  sulfate  and  is  important  for 
intracellular  mature  virion  morphogenesis  and  virus  infection  in  vitro 
and  in  vivo.  J  Virol  2000,  74:3353-3365. 

9.  Armstrong  JA,  Metz  DH,  Young  MR:  The  mode  of  entry  of  vaccinia  virus 
into  L  cells.  J  Gen  Virol  1973,  21:533-537. 

10.  Dales  S,  Kajioka  R:  the  cycle  of  multiplication  of  vaccinia  virus  in  earle's 
strain  I  cells.  I.  Uptake  and  penetration.  Virology  1964,  24:278-294. 

11.  Payne  LG,  Norrby  E:  Adsorption  and  penetration  of  enveloped  and  naked 
vaccinia  virus  particles.  J  Virol  1978,  27:19-27. 

12.  Vanderplasschen  A,  Hollinshead  M,  Smith  GL:  Intracellular  and  extracellular 
vaccinia  virions  enter  cells  by  different  mechanisms.  J  Gen  Virol  1998, 
79(Pt  4):877-887. 

13.  Munyon  W,  Paoletti  E,  Grace  JT  Jr:  RNA  polymerase  activity  in  purified 
infectious  vaccinia  virus.  Proc  Natl  Acad  Sci  USA  1967,  58:2280-2287. 

14.  Kates  JR,  McAuslan  BR:  Poxvirus  DNA-dependent  RNA  polymerase.  Proc 
Natl  Acad  Sci  USA  1967,  58:134-141. 

15.  Da  Fonseca  F,  Moss  B:  Poxvirus  DNA  topoisomerase  knockout  mutant 
exhibits  decreased  infectivity  associated  with  reduced  early 
transcription.  Proc  Natl  Acad  Sci  USA  2003,  100:1 1291-1 1296. 

16.  Werden  SJ,  Rahman  MM,  McFadden  G:  Poxvirus  host  range  genes.  Adv 
Virus  Res  2008,  71:135-171. 

17.  McFadden  G:  Poxvirus  tropism.  Nat  Rev  Microbiol  2005,  3:201-213. 

18.  Rubins  KH,  Hensley  LE,  Bell  GW,  Wang  C,  Lefkowitz  EJ,  Brown  PO, 

Reiman  DA:  Comparative  Analysis  of  Viral  Gene  Expression  Programs 
during  Poxvirus  Infection:  A  Transcriptional  Map  of  the  Vaccinia  and 
Monkeypox  Genomes.  PLoS  ONE  2008,  3:e2628. 

19.  Assarsson  E,  Greenbaum  JA,  Sundstrom  M,  Schaffer  L,  Hammond  JA, 
Pasquetto  V,  Oseroff  C,  Hendrickson  RC,  Lefkowitz  EJ,  Tscharke  DC,  et  al : 
Kinetic  analysis  of  a  complete  poxvirus  transcriptome  reveals  an 
immediate-early  class  of  genes.  Proc  Natl  Acad  Sci  USA  2008, 
105:2140-2145. 

20.  Guerra  S,  Lopez-Fernandez  LA,  Pascual-Montano  A,  Munoz  M,  Harshman  K, 
Esteban  M:  Cellular  gene  expression  survey  of  vaccinia  virus  infection  of 
human  HeLa  cells.  J  Virol  2003,  77:6493-6506. 

21.  Rubins  KH,  Hensley  LE,  Jahrling  PB,  Whitney  AR,  Geisbert  TW,  Huggins  JW, 
Owen  A,  Leduc  JW,  Brown  PO,  Reiman  DA:  The  host  response  to 
smallpox:  analysis  of  the  gene  expression  program  in  peripheral  blood 
cells  in  a  nonhuman  primate  model.  Proc  Natl  Acad  Sci  USA  2004, 
101:15190-15195. 

22.  Brum  LM,  Lopez  MC,  Varela  JC,  Baker  HV,  Moyer  RW:  Microarray  analysis  of 
A549  cells  infected  with  rabbitpox  virus  (RPV):  a  comparison  of  wild- 
type  RPV  and  RPV  deleted  for  the  host  range  gene,  SPI-1 .  Virology  2003, 
315:322-334. 

23.  Schuster  FL,  Visvesvara  GS:  Axenic  growth  and  drug  sensitivity  studies  of 
Balamuthia  mandrillaris,  an  agent  of  amebic  meningoencephalitis  in 
humans  and  other  animals.  J  Clin  Microbiol  1996,  34:385-388. 

24.  Tsuchiya  Y,  Tagaya  I:  Plaque  assay  of  variola  virus  in  a  cynomolgus 
monkey  kidney  cell  line.  Arch  Gesamte  Virusforsch  1970,  32:73-81. 

25.  Hui  EK,  Yap  EM,  An  DS,  Chen  IS,  Nayak  DP:  Inhibition  of  influenza  virus 
matrix  (Ml)  protein  expression  and  virus  replication  by  U6  promoter- 
driven  and  lentivirus-mediated  delivery  of  siRNA.  J  Gen  Virol  2004, 
85:1877-1884. 

26.  Lockhart  DJ,  Dong  H,  Byrne  MC,  Follettie  MT,  Gallo  MV,  Chee  MS, 

Mittmann  M,  Wang  C,  Kobayashi  M,  Horton  H,  Brown  EL:  Expression 
monitoring  by  hybridization  to  high-density  oligonucleotide  arrays.  Nat 

Biotechnol  1996,  14:1675-1680. 

27.  Willis  SN,  Chen  L,  Dewson  G,  Wei  A,  Naik  E,  Fletcher  Jl,  Adams  JM, 

Huang  DC:  Proapoptotic  Bak  is  sequestered  by  Mcl-1  and  Bcl-xL,  but  not 
Bcl-2,  until  displaced  by  BH3-only  proteins.  Genes  Dev  2005,  19:1294-1305. 

28.  Oda  E,  Ohki  R,  Murasawa  H,  Nemoto  J,  Shibue  T,  Yamashita  T,  Tokino  T, 
Taniguchi  T,  Tanaka  N:  Noxa,  a  BH3-only  member  of  the  Bcl-2  family  and 
candidate  mediator  of  p53-induced  apoptosis.  Science  2000, 
288:1053-1058. 

29.  Nakajima  W,  Tanaka  N:  Synergistic  induction  of  apoptosis  by  p53- 
inducible  Bcl-2  family  proteins  Noxa  and  Puma.  J  Nippon  Med  Sch  2007, 
74:148-157. 

30.  Debatin  KM,  Poncet  D,  Kroemer  G:  Chemotherapy:  targeting  the 
mitochondrial  cell  death  pathway.  Oncogene  2002,  21:8786-8803. 

31.  Pavlov  EV,  Priault  M,  Pietkiewicz  D,  Cheng  EH,  Antonsson  B,  Manon  S, 
Korsmeyer  SJ,  Mannella  CA,  Kinnally  KW:  A  novel,  high  conductance 


channel  of  mitochondria  linked  to  apoptosis  in  mammalian  cells  and 
Bax  expression  in  yeast.  J  Cell  Biol  2001,  155:725-731. 

32.  Shimizu  S,  Konishi  A,  Kodama  T,  Tsujimoto  Y:  BH4  domain  of  antiapoptotic 
Bcl-2  family  members  closes  voltage-dependent  anion  channel  and 
inhibits  apoptotic  mitochondrial  changes  and  cell  death.  Proc  Natl  Acad 
Sci  USA  2000,  97:3100-3105. 

33.  Bossy-Wetzel  E,  Green  DR:  Caspases  induce  cytochrome  c  release  from 
mitochondria  by  activating  cytosolic  factors.  J  Biol  Chem  1999, 
274:17484-17490. 

34.  Edick  MJ,  Tesfay  L,  Lamb  LE,  Knudsen  BS,  Miranti  CK:  Inhibition  of  integrin- 
mediated  crosstalk  with  epidermal  growth  factor  receptor/Erk  or  Src 
signaling  pathways  in  autophagic  prostate  epithelial  cells  induces 
caspase-independent  death.  Mol  Biol  Cell  2007,  18:2481-2490. 

35.  Melino  G,  Bernassola  F,  Ranalli  M,  Yee  K,  Zong  WX,  Corazzari  M,  Knight  RA, 
Green  DR,  Thompson  C,  Vousden  KH:  p73  Induces  apoptosis  via  PUMA 
transactivation  and  Bax  mitochondrial  translocation.  J  Biol  Chem  2004, 
279:8076-8083. 

36.  Earnshaw  WC,  Martins  LM,  Kaufmann  SH:  Mammalian  caspases:  structure, 
activation,  substrates,  and  functions  during  apoptosis.  Annu  Rev  Biochem 
1999,  68:383-424. 

37.  Orlowski  RZ:  The  role  of  the  ubiquitin-proteasome  pathway  in  apoptosis. 

Cell  Death  Differ  1999,  6:303-313. 

38.  Mailand  N,  Podtelejnikov  A V,  Groth  A,  Mann  M,  Bartek  J,  Lukas  J:  Regulation 
of  G(2)/M  events  by  Cdc25A  through  phosphorylation-dependent 
modulation  of  its  stability.  EMBO  J  2002,  21:591 1-5920. 

39.  Lindqvist  A,  Kallstrom  H,  Lundgren  A,  Barsoum  E,  Rosenthal  CK:  Cdc25B 
cooperates  with  Cdc25A  to  induce  mitosis  but  has  a  unique  role  in 
activating  cyclin  Bl-Cdkl  at  the  centrosome.  J  Cell  Biol  2005,  171:35-45. 

40.  Peng  CY,  Graves  PR,  Thoma  RS,  Wu  Z,  Shaw  AS,  Piwnica-Worms  H:  Mitotic 
and  G2  checkpoint  control:  regulation  of  14-3-3  protein  binding  by 
phosphorylation  of  Cdc25C  on  serine-216.  Science  1997,  277:1501-1505. 

41.  Bowe  DB,  Kenney  NJ,  Adereth  Y,  Maroulakou  IG:  Suppression  of  Neu- 
induced  mammary  tumor  growth  in  cyclin  D1  deficient  mice  is 
compensated  for  by  cyclin  E.  Oncogene  2002,  21:291-298. 

42.  Bales  E,  Mills  L,  Milam  N,  McGahren-Murray  M,  Bandyopadhyay  D,  Chen  D, 
Reed  JA,  Timchenko  N,  van  den  Oord  JJ,  Bar-Eli  M,  et  at.  The  low  molecular 
weight  cyclin  E  isoforms  augment  angiogenesis  and  metastasis  of 
human  melanoma  cells  in  vivo.  Cancer  Res  2005,  65:692-697. 

43.  Deane  NG,  Lee  H,  Hamaamen  J,  Ruley  A,  Washington  MK,  LaFleur  B, 
Thorgeirsson  SS,  Price  R,  Beauchamp  RD:  Enhanced  tumor  formation  in 
cyclin  D1  x  transforming  growth  factor  betal  double  transgenic  mice 
with  characterization  by  magnetic  resonance  imaging.  Cancer  Res  2004, 
64:1315-1322. 

44.  Nakamoto  M:  Eph  receptors  and  ephrins.  Int  J  Biochem  Cell  Biol  2000, 
32:7-12. 

45.  Ewing  RM,  Chu  P,  Elisma  F,  Li  H,  Taylor  P,  Climie  S,  McBroom-Cerajewski  L, 
Robinson  MD,  O'Connor  L,  Li  M,  et  al:  Large-scale  mapping  of  human 
protein-protein  interactions  by  mass  spectrometry.  Mol  Syst  Biol  2007, 
3:89. 

46.  Latella  L,  Sacco  A,  Pajalunga  D,  Tiainen  M,  Macera  D,  D'Angelo  M,  Felici  A, 
Sacchi  A,  Crescenzi  M:  Reconstitution  of  cyclin  D1 -associated  kinase 
activity  drives  terminally  differentiated  cells  into  the  cell  cycle.  Mol  Cell 
Biol  2001,  21:5631-5643. 

47.  Siegert  JL,  Rushton  JJ,  Sellers  WR,  Kaelin  WG  Jr,  Robbins  PD:  Cyclin  D1 
suppresses  retinoblastoma  protein-mediated  inhibition  of  TAFII250 
kinase  activity.  Oncogene  2000,  19:5703-5711. 

48.  King  KL,  Cidlowski  JA:  Cell  cycle  regulation  and  apoptosis.  Annu  Rev 
Physiol  1998,  60:601-617. 

49.  Martin  LG,  Demers  GW,  Galloway  DA:  Disruption  of  the  G1/S  transition  in 
human  papillomavirus  type  16  E7-expressing  human  cells  is  associated 
with  altered  regulation  of  cyclin  E.  J  Virol  1998,  72:975-985. 

50.  Hallstrom  TC,  Nevins  JR:  Balancing  the  decision  of  cell  proliferation  and 
cell  fate.  Cell  Cycle  2009,  8:532-535. 

51.  Roussel  MF:  The  INK4  family  of  cell  cycle  inhibitors  in  cancer.  Oncogene 
1999,  18:5311-5317. 

52.  Kohn  KW:  Molecular  interaction  map  of  the  mammalian  cell  cycle 
control  and  DNA  repair  systems.  Mol  Biol  Cell  1999,  10:2703-2734. 

53.  Falck  J,  Mailand  N,  Syljuasen  RG,  Bartek  J,  Lukas  J:  The  ATM-Chk2-Cdc25A 
checkpoint  pathway  guards  against  radioresistant  DNA  synthesis.  Nature 
2001,  410:842-847. 


Alkhalil  et  al.  Virology  Journal  2010,  7:173 
http://www.virologyj.eom/content/7/1/1 73 


Page  19  of  19 


54.  Mailand  N,  Falck  J,  Lukas  C,  Syljuasen  RG,  Welcker  M,  Bartek  J,  Lukas  J: 

Rapid  destruction  of  human  Cdc25A  in  response  to  DNA  damage. 

Science  2000,  288:1425-1429. 

55.  Lew  J:  MAP  kinases  and  CDKs:  kinetic  basis  for  catalytic  activation. 

Biochemistry  2003,  42:849-856. 

56.  Sherr  G,  Roberts  JM:  CDK  inhibitors:  positive  and  negative  regulators  of 
Gl-phase  progression.  Genes  Dev  1999,  13:1501-1512. 

57.  Donovan  JC,  Rothenstein  JM,  Slingerland  JM:  Non-malignant  and  tumor- 
derived  cells  differ  in  their  requirement  for  p27Kip1  in  transforming 
growth  factor-beta-mediated  G1  arrest.  J  Biol  Chem  2002, 

277:41686-41692. 

58.  Manoukian  AS,  Woodgett  JR:  Role  of  glycogen  synthase  kinase-3  in 
cancer:  regulation  by  Wnts  and  other  signaling  pathways.  Adv  Cancer  Res 
2002,  84:203-229. 

59.  Sakamoto  K,  Creamer  BA,  Triplett  AA,  Wagner  KU:  The  Janus  kinase  2  is 
required  for  expression  and  nuclear  accumulation  of  cyclin  D1  in 
proliferating  mammary  epithelial  cells.  Mol  Endocrinol  2007,  21:1877-1892. 

60.  Goldstone  S,  Pavey  S,  Forrest  A,  Sinnamon  J,  Gabrielli  B:  Cdc25-dependent 
activation  of  cyclin  A/cdk2  is  blocked  in  G2  phase  arrested  cells 
independently  of  ATM/ATR.  Oncogene  2001,  20:921-932. 

61.  Kumagai  A,  Dunphy  WG:  The  cdc25  protein  controls  tyrosine 
dephosphorylation  of  the  cdc2  protein  in  a  cell-free  system.  Cell  1991, 
64:903-914. 

62.  Lavin  MF,  Shiloh  Y:  The  genetic  defect  in  ataxia-telangiectasia.  Annu  Rev 
Immunol  1997,  15:177-202. 

63.  Canman  CE,  Lim  DS,  Cimprich  KA,  Taya  Y,  Tamai  K,  Sakaguchi  K,  Appella  E, 
Kastan  MB,  Siliciano  JD:  Activation  of  the  ATM  kinase  by  ionizing 
radiation  and  phosphorylation  of  p53.  Science  1998,  281:1677-1679. 

64.  Cliby  WA,  Roberts  CJ,  Cimprich  KA,  Stringer  CM,  Lamb  JR,  Schreiber  SL, 
Friend  SH:  Overexpression  of  a  kinase-inactive  ATR  protein  causes 
sensitivity  to  DNA-damaging  agents  and  defects  in  cell  cycle 
checkpoints.  EMBO  J  1998,  17:159-169. 

65.  Zhao  H,  Piwnica-Worms  H:  ATR-mediated  checkpoint  pathways  regulate 
phosphorylation  and  activation  of  human  Chkl.  Mol  Cell  Biol  2001, 
21:4129-4139. 

66.  Bahassi  el  M,  Conn  CW,  Myer  DL,  Flennigan  RF,  McGowan  CH,  Sanchez  Y, 
Stambrook  PJ:  Mammalian  Polo-like  kinase  3  (Plk3)  is  a  multifunctional 
protein  involved  in  stress  response  pathways.  Oncogene  2002, 
21:6633-6640. 

67.  Sanchez  Y,  Wong  C,  Thoma  RS,  Richman  R,  Wu  Z,  Piwnica-Worms  H, 

Elledge  SJ:  Conservation  of  the  Chkl  checkpoint  pathway  in  mammals: 
linkage  of  DNA  damage  to  Cdk  regulation  through  Cdc25.  Science  1997, 
277:1497-1501. 

68.  Hengartner  MO:  The  biochemistry  of  apoptosis.  Nature  2000,  407:770-776. 

69.  Everett  H,  McFadden  G:  Poxviruses  and  apoptosis:  a  time  to  die.  Curr  Opin 
Microbiol  2002,  5:395-402. 

70.  Thome  M,  Schneider  P,  Hofmann  K,  Fickenscher  H,  Meinl  E,  Neipel  F, 
Mattmann  C,  Burns  K,  Bodmer  JL,  Schroter  M,  et  al :  Viral  FLICE-inhibitory 
proteins  (FLIPs)  prevent  apoptosis  induced  by  death  receptors.  Nature 
1997,  386:517-521. 

71.  Shisler  JL,  Senkevich  TG,  Berry  MJ,  Moss  B:  Ultraviolet-induced  cell  death 
blocked  by  a  selenoprotein  from  a  human  dermatotropic  poxvirus. 

Science  1998,  279:102-105. 

72.  Komiyama  T,  Ray  CA,  Pickup  DJ,  Howard  AD,  Thornberry  NA,  Peterson  EP, 
Salvesen  G:  Inhibition  of  interleukin-1  beta  converting  enzyme  by  the 
cowpox  virus  serpin  CrmA.  An  example  of  cross-class  inhibition.  J  Biol 
Chem  1994,  269:19331-19337. 

73.  Quan  LT,  Caputo  A,  Bleackley  RC,  Pickup  DJ,  Salvesen  GS:  Granzyme  B  is 
inhibited  by  the  cowpox  virus  serpin  cytokine  response  modifier  A.  J  Biol 
Chem  1995,  270:10377-10379. 

74.  Dobbelstein  M,  ShenkT:  Protection  against  apoptosis  by  the  vaccinia 
virus  SPI-2  (B13R)  gene  product.  J  Virol  1996,  70:6479-6485. 

75.  Davies  MV,  Chang  HW,  Jacobs  BL,  Kaufman  RJ:  The  E3L  and  K3L  vaccinia 
virus  gene  products  stimulate  translation  through  inhibition  of  the 
double-stranded  RNA-dependent  protein  kinase  by  different 
mechanisms.  J  Virol  1993,  67:1688-1692. 

76.  Wasilenko  ST,  Stewart  TL,  Meyers  AF,  Barry  M:  Vaccinia  virus  encodes  a 
previously  uncharacterized  mitochondrial-associated  inhibitor  of 
apoptosis.  Proc  Natl  Acad  Sci  USA  2003,  100:14345-14350. 


77.  Parrino  J,  Hotchkiss  RS,  Bray  M:  Prevention  of  immune  cell  apoptosis  as 
potential  therapeutic  strategy  for  severe  infections.  Emerg  Infect  Dis  2007, 
13:191-198. 

78.  Afonso  CL,  Tulman  ER,  Lu  Z,  Zsak  L,  Kutish  GF,  Rock  DL:  The  genome  of 
fowlpox  virus.  J  Virol  2000,  74:3815-3831. 

79.  Soloski  MJ,  Holowczak  JA:  Characterization  of  supercoiled  nucleoprotein 
complexes  released  from  detergent-treated  vaccinia  virions.  J  Virol  1981, 
37:770-783. 

80.  Puvion-Dutilleul  F,  Leduc  EH,  Puvion  E:  The  osmium  ammine-S02  staining 
method  for  studying  the  in  situ  configuration  of  viral  genomes  in 
ultrathin  sections  of  DNA  virus  infected  cells.  Biol  Cell  1996,  87:133-141. 

81.  Lilley  CE,  Chaurushiya  MS,  Weitzman  MD:  Chromatin  at  the  intersection  of 
viral  infection  and  DNA  damage.  Biochim  Biophys  Acta  2009,  1799:319-327. 

82.  Fanning  E,  Zhao  K:  SV40  DNA  replication:  from  the  A  gene  to  a 
nanomachine.  Virology  2009,  384:352-359. 

83.  Scaife  RM,  Job  D,  Langdon  WY:  Rapid  microtubule-dependent  induction 
of  neurite-like  extensions  in  NIH  3T3  fibroblasts  by  inhibition  of  ROCK 
and  Cbl.  Mol  Biol  Cell  2003,  14:4605-4617. 

84.  McGlade  J,  Brunkhorst  B,  Anderson  D,  Mbamalu  G,  Settleman  J,  Dedhar  S, 
Rozakis-Adcock  M,  Chen  LB,  Pawson  T:  The  N-terminal  region  of  GAP 
regulates  cytoskeletal  structure  and  cell  adhesion.  EMBO  J  1993, 
12:3073-3081. 

85.  Pawlak  G,  Helfman  DM:  Post-transcriptional  down-regulation  of  ROCKI/ 
Rho-kinase  through  an  MEK-dependent  pathway  leads  to  cytoskeleton 
disruption  in  Ras-transformed  fibroblasts.  Mol  Biol  Cell  2002,  13:336-347. 

86.  Giancotti  FG,  Ruoslahti  E:  Integrin  signaling.  Science  1999,  285:1028-1032. 

87.  Pochynyuk  0,  Stockand  JD,  Staruschenko  A:  Ion  channel  regulation  by  Ras, 
Rho,  and  Rab  small  GTPases.  Exp  Biol  Med  (Maywood)  2007,  232:1258-1265. 

88.  Everett  H,  Barry  M,  Sun  X,  Lee  SF,  Frantz  C,  Berthiaume  LG,  McFadden  G, 
Bleackley  RC:  The  myxoma  poxvirus  protein,  M11L,  prevents  apoptosis 
by  direct  interaction  with  the  mitochondrial  permeability  transition 
pore.  J  Exp  Med  2002,  196:1 127-1 139. 

89.  Schubert  U,  Ferrer-Montiel  AV,  Oblatt-Montal  M,  Henklein  P,  Strebel  K, 
Montal  M:  Identification  of  an  ion  channel  activity  of  the  Vpu 
transmembrane  domain  and  its  involvement  in  the  regulation  of  virus 
release  from  H I V-1 -infected  cells.  FEBS  Lett  1996,  398:12-18. 

90.  Desai  SA,  Bezrukov  SM,  Zimmerberg  J:  A  voltage-dependent  channel 
involved  in  nutrient  uptake  by  red  blood  cells  infected  with  the  malaria 
parasite.  Nature  2000,  406:1001-1005. 

91.  Hrimech  M,  Yao  XJ,  Branton  PE,  Cohen  EA:  Human  immunodeficiency 
virus  type  1  Vpr-mediated  G(2)  cell  cycle  arrest:  Vpr  interferes  with  cell 
cycle  signaling  cascades  by  interacting  with  the  B  subunit  of  serine/ 
threonine  protein  phosphatase  2A.  EMBO  J  2000,  19:3956-3967. 

92.  He  J,  Choe  S,  Walker  R,  Di  Marzio  P,  Morgan  DO,  Landau  NR:  Human 
immunodeficiency  virus  type  1  viral  protein  R  (Vpr)  arrests  cells  in  the 
G2  phase  of  the  cell  cycle  by  inhibiting  p34cdc2  activity.  J  Virol  1995, 
69:6705-6711. 

93.  Re  F,  Braaten  D,  Franke  EK,  Luban  J:  Human  immunodeficiency  virus  type 
1  Vpr  arrests  the  cell  cycle  in  G2  by  inhibiting  the  activation  of 
p34cdc2-cyclin  B.  J  Virol  1995,  69:6859-6864. 

94.  Zhao  RY,  Elder  RT:  Viral  infections  and  cell  cycle  G2/M  regulation.  Cell  Res 
2005,  15:143-149. 

95.  Janssens  V,  Goris  J:  Protein  phosphatase  2A:  a  highly  regulated  family  of 
serine/threonine  phosphatases  implicated  in  cell  growth  and  signalling. 

Biochem  J  2001,  353:417-439. 


doi:1 0.1 1 86/1 743-422X-7-1 73 

Cite  this  article  as:  Alkhalil  et  al.:  Gene  expression  profiling  of 
monkeypox  virus-infected  cells  reveals  novel  interfaces  for  host-virus 
interactions.  Virology  Journal  2010  7:173. 


